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R-685 


PRELIMINARY SOFTWARE DEVELOPMENT 
FOR OPTICAL MIRROR FIGURE CONTROL 


ABSTRACT 

The maintenance of accurate primary mirror figure in the face of 
environmental disturbances is the key to the achievement of diffraction- 
limited performance in a large space telescope. In order to develop the 
concepts of optical mirror figure control, an experimental program has 
been initiated at the Marshall Space Flight Center, Huntsville, Alabama. 

A major component in this experiment will be an XDS Sigma 5/7 digital 
computer which will realize the control algorithm. This report presents 
a description of a software package which realizes linear optimal, 
simplified linear and iterative optimal control algorithms. The software, 
in addition, provides for interactive communication between the operator 
and the computer, and interaction between the computer and the exper- 
imental hardware elements. A brief description of a small hybrid computer 
system is also presented. 
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CHAPTER 1 

AN EXPERIMENTAL ACTIVE MIRROR 


1.0 Introduction 


Astronomical observations through a large earth-based telescope suffer 
from limitations placed on the resolving power of the telescope by fluctua- 
tions in the earth' s atmosphere. As part of a space orbiting astronomical 
laboratory, however, a telescope would not be subject to these limitations. 

A large instrument which is diffraction - limited over a major part of its 
useful spectrum of observation was envisioned. Maximum resolving 

power requires extremely accurate maintenance of the figure of the 
1 2 

primary mirror. • 

Although it was possible to suitably polish such a large mirror to the 

desired surface accuracy, stresses introduced by thermal variations 

in the mirror and fluctuations in support structure loads, structural 

instability and the elimination of gravity loading in orbit could create 

surface perturbations which would exceed the surface accuracy limits 

required for diffraction limited performance. As a result, investigators 

have attempted to develop techniques for actively correcting the mirror 

figure in a space environment and a number of promising control tech- 

1 _ 8 

niques have been developed. The development and application of 
these control concepts is one of the key challenges facing the designer 
of the large space telescope. 

While the large scale digital computer is an extremely valuable tool 
for the analysis and simulation of complicated mirror figure control 
systems, the results obtained are only as reliable as the modelling 
accuracy of the physical components in the system. Accurate modelling 
requires a considerable amount of intuition if the trade-off between modelling 
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accuracy and computation time is to be satisfactorily resolved. 

Often, terms neglected in the modelling process are of key importance 
to the overall system design. 


To resolve these problems it is important to have some way of checking 
the results of numerical analysis against actual system behavior. Such 
checks are furnished by an experimental program. 

Experimental work in the past has been largely conducted using analog 
devices to synthesize figure actuator commands from surface error 
measurements. As a result of the expense and time associated with 
programming a general purpose analog computer or constructing a 
special purpose analog system, it has been difficult to explore the full 
spectrum of control solutions or efficiently process experimental data. 

In response to these problems, a more efficient experimental tool 
has evolved in the hybrid digital analog computer system. Spurred 
by declining cost hybrid computing systems are appearing in a wide 
variety of laboratory environments. This report describes the application 
of a hybrid computing system to an experimental active mirror and the 
required software packages, 

1. 1 Experiment Design Concept 


The experimental system consists of the primary mirror fitted 
with actuators for figure modification, a mirror figure sensor and 
an "on line" control system processing the figure errors measured 
by the sensor to provide proper corrective signals to the actuators. 


2 



Following a modern approach a digital computer has been selected 
as the control processor of the experimental system. The utilization 
of a general purpose computer has a number of advantages, including: 

1. Programmability permitting a large number of different 
primary mirror control configurations to be investigated 
without extensive hardware modification. 

2. The ability to handle a number of auxilliary tasks such as 
experimental data processing and display. 

3. Characteristics similar to the system computer which will 
be used in an orbiting astronomical observatory. 

1. 2 A Computer Controlled Active Mirror 

A block diagram of an experimental active mirror is shown in Figure 1,2, 1. 
The digital computer consists of a central processor, a random access core 
memory and an input-output processor. The central processor handles arith- 
metic operations, logic operations and some date transfer operations using 
instructions extracted from the core memory. The input-output processor 
controls the transfer of information from the central processor and core 
memory to computer peripherals which are part of the interface with 
the real world. The core memory holds two types of stored information, 
program instructions and program storage. The program instructions 
tell the computer what to do with information extracted from program 
storage and input devices. 

The interface consists of a number of components. The digital output 
buffers are a set of addressable registers which temporarily store digital 
information which is used by devices such as the external switch 
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and data displays. The buffer registers temporarily hold digital data 
which is converted to an analog signal by the digital to analog converters. 
The output signals from the digital to analog converters are inputs to 
display devices such as meters and scopes and an external switch 
which expands the capacity of the converter by permitting a larger number 
of channels to be addressed at a lower switching rate. The external 
switch is controlled by digital signals from the digital output buffer. 

Analog signals to be entered into the computer pass through a computer- 
controlled switch or multiplexer; the single output of the multiplexer 
is converted to a digital signal by the analog to digital converter. 

The digital signal is then transferred to the central processor or the 
correct location in the core memory by the input-output processor. 

External to the computer and its interface is the mirror figure sensor . 
Signals from the figure sensor pass through an analog signal processor , 
if required, to the multiplexer. The outputs from the external switch 
pass to a power amplifier which supplies the energy required to drive the 
mirror figure actuators. 
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CHAPTER 2 

GENERAL FEATURES OF THE EAM SOFTWARE PACKAGE 


2. 0 General Specifications 

The software package is divided into 4 major sections: 


1. EAM Package for Linear Optimal or Simplified Linear 
Figure Control (EAM). 

2. Parameter Optimization Package for Iterative 
Figure Control (POP). 

3. Mathematical Operations Package (MOP). 

4. Input Output Operations Package (IOP). 


The EAM package contains programs which realize the Linear Optimal 
(LOCS) and Simplified Linear (SLCS) control laws. The programs also 
calculate system gains and initialize system parameters. An operating 
sequence monitoring capability is also provided which prevents operator 
instruction sequence errors. The programs also permit the operator 
to operate and diagnose the EAM remotely using a typewriter console. 
The subroutines also realize servo loops for the actuators. The addition 
of actuator output control loops reduces the sensitivity of the control 
system to variations in actuator scale factors due to nonlinearity and 
structural instability. 


The Parameter Optimization Subroutine Package provides the software 
for the Iterative Figure Control System (IFCS). The gradient of the 
performance index 


/ 

J = Xj 

where x^ is the measured figure error vector 


( 2 . 0 . 1 ) 
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is generated by perturbing the elements of the actuator servo input 
vector. The resultant gradient information is used in a steepest descent 
procedure to generate a sequence of actuator outputs which approach 
the linear optimal control solution. Modified versions of the steepest 
descent algorighm featuring improved rates of convergence may also 
be selected by the operator. 

The Mathematical Subroutine Package is used by all the EAM programs 
to perform mathematical operations such as matrix addition, subtraction, 
multiplication and inversion, etc. All matrices are stored as single 
dimensioned arrays. While this complicates the programming somewhat 
the problem of carrying dimensions through subroutines and restrictions 
on variable dimensioning are eliminated. 

The Input-Output Subroutine Package has been developed at the Charles 
Stark Draper Laboratory to permit highly efficient handling of input 
data and output print in terms of core storage for program instructions. 
Further savings are realized by reading and printing headings associated 
with input data. This procedure eliminates many format statements 
which would otherwise require a large amount of core. 

2. 1 Computer-Operator Communication 


The computer exists in an isolated electronic world. All communication 
with the external environment must take place through channels which 
process the information in a form suitable for computer or peripheral 
use. The data channels include analog to digital converters and digital 
information buffers for input and digital to analog converters and analog 

buffers for output. Two levels of software are envisioned for effective 
communication. 
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The first division is a machine oriented set or programs which perform 
data transfer on the lowest channel level. Such software is written in 
machine language and provides an interface between the higher level 
FORTRAN programs and computer hardware. For example, if it 
was desired to plot information on an analog X-Y plotter the computer 
programmer would call a FORTRAN program PLOT which would scale 
the data and provide interpolated values if necessary. The resultant 
scaled data would be transfered to the analog channels assigned to the 
X and Y coordinates by a machine language subroutine which would 
receive the channel numbers and the data from the FORTRAN program. 
Such a procedure offers a maximum amount of flexibility since PLOT 
in FORTRAN is easily modified or expanded while the machine oriented 
subroutine can remain fixed. In addition, conversion to a new computer 
facility can generally be accomplished with minor changes in the FORTRAN 
routines and the substitution of a new set of machine language subroutines. 

The EAM software is designed to take maximum advantage of the charac- 
teristics of a variety of available computer peripherals. Typical input- 
output peripherals may be classified as shown in table 2. 1-1. The 
line printer, high speed card reader and magnetic tape are ideal devices 
for transferring large blocks of data. A typewriter or teletype, on the 
other hand, is best suited for making modifications to or extracting 
small blocks of data from the program data base. 

2. 2 The Data Buss Concept 


In order to assure maximum access to program storage a "data buss" 
has been established. The "data buss" is a one dimensional array 
designated XV. All program variables are equivalenced to sections 
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of XV. XV is transmitted throughout the EAM software by means of 
a labelled common block BLKMFC. Transmission of information via 
COMMON storage eliminates the requirement for long parameter lists 
in subroutine calls. 


Table 2. 1-1 Peripheral device characteristics. 


DEVICE DATA TRANSFER RATE COST SIZE 


LOW 

HIGH 

LOW 

HIGH 

SMALL 

LARGE 


/ 


/ 


/ 


/ 


/ 


/ 
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/ 
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CHAPTER 3 


DESCRIPTION OF PRINCIPAL EAM SOFTWARE COMPONENTS 


3. 0 Introduction 


The following sections of Chapter 3 are general descriptions of the 
major components of the EAM software package. User descriptions 
are presented in Chapter 4 while program listings are contained in 
Appendices A, B, C, and D. The software is written using a combi- 
nation of XDS Fortran II and IV and is designed for execution on an 
XDS Sigma 5/7 computer.^ 11 

Program and subprogram names (MAIN, TYPCON, ACTCAL, etc. ) 
and program variables (MODOP, MODE, XV, etc. ) are capitalized in 
the text. Indexed variable names are usually followed by the index in 
brackets (XV(I), UFV(I), etc. ). 

3. 1 Main Program for the Experimental Active Mirror 


Program MAIN provides the skeletal structure for the EAM. The 
program initially interrogates the typewriter through TYPCON. Control 
remains with the typewriter until the control system operating mode 
(LOCS, SLCS, or IFCS) is defined. At that point the defined operating 
mode (MODOP) is stored and MAIN monitors the sequence of typewriter 
commands to the system. The correct operating sequence is: 


1. Identify control configuration (MODE = MODOP > 5). 

2. Initialize EAM (MODE = 1). 

3. Start EAM (MODE = 2). 


Operation of the EAM control system may be terminated at any time by 
setting sense switch NSNSWT (defined by input data) which transfers 
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control back to the typewriter. The typewriter may then be used to 
examine or change any piece of information in the DATA BUSS. When 
control is returned to MAIN the status of MODE and MODOP is checked 
to ensure that they have not been altered. Alteration generates a 
typed operator error message and returns control to the typewriter for 
correction. A flow diagram of the main program is shown in Fig. 3. 1. 1. 

3. 2 Typewriter Control of the EAM 

Since the computer is generally some distance from the experiment 
it is desirable to provide a remote control capability adjacent to the 
experiment. In the current configuration this capability is provided by 
a typewriter which communicates with the EAM software through the 
subroutine TYPCON. 

When TYPCON is called with NENTRY = 2 the typewriter is interrogated 
for a value of MODE. MODE can have integer values between 1 and 10. 
The operations corresponding to these modes are 

MODE = 1 Initialize the Mirror Figure Control System (MFCS) 
Control is transferred to MAIN which reads in data 
and initializes the EAM providing MODOP has been 
defined. If MODOP was not defined an error message 
is typed and new value of MODE requested by TYPCON. 

MODE = 2 Start the Mirror Figure Control System. Control is 
transferred to MAIN which starts the control system 
providing the previous value of MODE was 1 and 
MODOP has not changed since MODE was equal to 1. 

If these conditions are not satisfied an error message 
is typed and a new MODE requested. 
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MODE = 3 


MODE = 4 


MODE = 5 


MODE =6 


Diagnostic mode. The abbreviated name of a variable 
in the DATA BUSS is requested. The typed reply is 
tested for conformity with the cataloged names. If 
the reply is not in the catalog an error message is 
typed and a new name requested. If the name is accep- 
ted, it is typed. The subscript(s) of subscripted variables 
are then requested. The operator types in the subscript(s) 
which is retyped by the computer. TYPCON then inter- 
rogates the value of the identified parameter and types 
it. A new variable name is then requested. Diagnosis 
may be continued at this point by typing in a variable 
name or terminated by typing DEND. If DEND is typed 
TYPCON will request a new value for MODE. To save 
time and storage space the variable names are abbreviated 
as indicated in Table 3. 2. 1. 

Test actuators . Subroutine ACTCAL is called and the 
actuator scale factors printed on the line printer. 

The end of ACTCAL is signaled by a typed message 
and a new value for MODE is requested. 

Test mirror. Subroutine MIRCAL is called and the 
reduced deformation - actuator input matrix generated. 

The resultant matrix is printed and the end of calibra- 
tion indicated by a typed message and the request 
for a new MODE. 

Simplified linear figure control. TYPCON sets 

MODOP = 6 and transfers control to MAIN. This 

operation defines a simplified linear control mode 

which may then be initialized and operated. 
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Fig. 3. 1. 1 Main program for the experimental 
active mirror (MAIN). 
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* TYPEWRITER INTERROGATION 
FOR CONTROL INSTRUCTION 


. 3. 1. 1 Main program for the experimental 
active mirror (MAIN). (Cont) 
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Table 3. 2. 1 Program variable abbreviations 


PROGRAM VARIABLE 

ABBREVIATED NAME 

AIM 

AIM 

AM 

AM 

ASCALV 

ASCV 

DUMV 

DUMV 

FSCALV 

FSCV 

GAINV 

GANV 

GAINM 

GANM 

LACTV 

LACV 

UFA V 

UFAV 

UFV 

UFV 

XFRV 

XFR V * 

XFSV 

XFSV 

XFV 

XFV 

XV 

XV 

YFSV 

YFSV 
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MODE = 7 


MODE = 8 


MODE = 9 


Linear optimal figure control. TYPCON sets MODOP 
and transfers control to MAIN. MODOP = 7 specifies 
a linear optimal control mode which may then be ini- 
tialized and operated. 

Iterative figure control . MODOP is set equal to 8 by 
TYPCON. MODOP = 8 defines an iterative control 
mode which may then be initialized and operated. 

Modify data buss value. This mode establishes a 
condition in which it is possible for the operator to 
modify any element in the data buss. The element is 
identified and indexed using the procedure under 
MODE = 3. A new value for the completely identified 
variable is then accepted; the variable value changed 
and the resulting value printed. 


A flow diagram of TYPCON is shown in Fig. 3. 2. 1. 



ENTER 



Fig. 3. 2. 1 Flow diagram of the typewriter control algorithm 
for the experimental active mirror (TYPCON). 
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Fig. 3. 2.1 Flow diagram of the typewriter control algorithm for 
the experimental active mirror (TYPCON). (Cont) 
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3 . 3 Experimental Active Mirror Control Algorithm 

The mirror figure control algorithm is realized in subroutine MFCS. 

When the parameter NENTRY = 1 the gains associated with the control 
law and the mirror actuators are read in and the feedback gain matrices 
for the Linear Optimal Control System or Simplified Linear Control 
System calculated. If Iterative Figure Control is desired the main 
parameter optimization program MINFCN is initialized. If NENTRY = 2 
the program calculates the linear optimal or simplified linear actuator 
commands. If NENTRY = 4 MFCS transfers the commands calculated by 
the Iterative Figure Control System to the actuator servos. A flow diagram 
of MFCS is shown in Fig. 3.3.1. The computational flow of the Iterative 
Figure Control System (MINFCN) is shown in Fig. 3.3.2. A detailed 
discription of the algorithm is presented in Ref. 1. 

3 . 4 Experimental Active Mirror Initialization 

The basic initialization of the EAM is accomplished by subroutine INIT. 

The initialization operations depend on the value of the parameter NENTRY 
which is determined when INIT is called. If NENTRY is equal to 1 INIT 
reads in the basic parameters associated with the EAM hardware components. 
When NENTRY is equal to 2 INIT initializes the program variables such 
as figure actuator output, input, and the figure measurement vector 
associated with the figure sensor. 

3 . 5 Actuator Command Processing 


Subroutine ACTCMD is called with an integer parameter I which identi- 
fies the Ith actuator. The subprogram interrogates the Ith actuator 
output sensor through machine language coding and transfers the reading 
to the calling program using UFAV(I). The Ith actuator input command 
UFV(I) is then multiplied by a scale factor ASCALV(I) which permits 
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ENTER 


RETURN 



RETURN 


Fig. 3.3.1 Flow diagram of mirror figure control 
system subroutine (MFCS). 
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compensation of scale differences between different actuators. The 
scaled actuator command is transferred to the actuator power modulator 
by a set of machine language statements which communicate with the 
D/A or D/D actuator channels. The subroutine features built-in adjus- 
table time delays which allow sufficient time for the interface components 
to operate. 

3. 6 Figure Sensor Measurement Processing 

Figure error measurements are obtained by the computer using subroutine 
FIGSEN. FIGSEN supplies the X and Y coordinates (XFSV(I), YFSV(I)) 
of the I' th figure error XFV(I) to the figure sensor image dissector. 

The figure sensor error measurement is returned to FIGSEN and mul- 
tiplied by the scale factor FSCALV<I). In order to reduce the effects 
of sensor noise the measurements are repeated NMEAS times and the 
results averaged. NMEAS is a local variable in FIGSEN defined when 
FIGSEN is called with NENTRY = 1. If the noise components associated 
with subsequent measurements are uncorrelated the standard deviation 
of the measurement is reduced (NMEAS)y times by the averaging process. 

3. 7 Actuator Calibration 


The scale factors associated with the I' th mirror figure actuator is de- 
termined by perturbing the actuator input UFV(I) by an amount DACT and 
observing the actuator output UFAV(I) measured by a sensor. The scale 
DUMV(I) is calculated by the expression 


(UFAV(I)) ufv(i) = DACT " (UFAV(I)) UFV(I) = - DACT 
DXJMV(I) = — 2 DACT 


( 3. 7. 1) 
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In order to minimize the effect of noise the measurements and scale 
computation are repeated NMEAS times and the results averaged. 

NMEAS is read in together with DACT. The subroutine prints the vector 
DUMV when all the scales factors have been determined. 

3. 8 Mirror Calibration 


The reduced deformation - force matrix associated with the mirror is 
generated by perturbing the actuator outputs and measuring the resulting 
figure changes. The actuator command UFV(J) is modified and the resulting 
figure errors XFV(I) I = 1, N observed. The sensitivity coefficients are 
calculated from the equation 


AM(I, J) 


(XFV(I)) UFV(J) ■ DACT - (XFV(I)) U.FV(J) 

2 DACT 


- DACT 


(3. 8-1) 


where AM(I, J) is the I, J th element of the deformation force matrix AM. 


In order to minimize the effects of figure sensor noise the test is repeated 
NMEAS times and the results averaged. NMEAS and DACT are read 
in by MIRCAL. If the noise components in subsequent measurements 
are uncorrelated the averaging process reduces the measurement 
deviation (NMEAS )y* times. 
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CHAPTER 4 

USER DESCRIPTIONS OF EAM SOFTWARE COMPONENTS 


4„ 0 Introduction 


The following pages contain user descriptions of the software components 
of the EAM. The descriptions list: 

1. Program name. 

2. The purpose of the program. 

3. The calling statement form. 

4. Definitions of the program parameters. 

5. Card input data description. 

The descriptions are divided into 4 sections: 

1. EAM Package (EAM). 

2. Parameter Optimization Package (POP). 

3. Mathematical Operations Package (MOP). 

4. Input Output Operations Package (IOP). 

Listings of the programs in each section are contained in appendices 
A, B, C and D respectively. 

The notation RBR in the data description indicates that the data is 
for a doubly dimensional array and is read in row by row. 

Descriptions are provided for the following subprograms: 
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EXPERIMENTAL ACTIVE MIRROR 
SOFTWARE COMPONENTS 


EAM PACKAGE 

MAIM MAIN PROGRAM FOR THF EXPFR I MENTAL ACTIVE MIRROR 

ACTCAl. ACTUATOR CALIBRATION 

ACIC.MO SCALES AND TRANSFERS ACTUATOR COMMANDS AND ACTUATOR OUTPUT 

MEASUREMENTS 

CHNIG MODIFIES THE VALUE OF A VARIABLE 

FIGSFN MEASURES THE FIGURE ERROR AT A DISCRETE POINT ON THE REFLECTING 

SURFACE OF THE MIRROR 

INIT EXPERIMENTAL ACTIVE MIRROR INITIALIZATION 

MFCS REALIZFS MIRROR FIGURE CONTROL SYSTEMS 

MIRCAL MIRROR CALIBRATION 

MODCHK CHECKS VALUES OF MODE AND MOOOP 

REALT INTERROGATES REAL TIME CLOCK 

RFDUAM GENERATES AR AMD ARR 

TYPE. ON TYPEWRITER CONTROL 


PARAMETER OPTIMIZATION PACKAGE 

FPCHNG STEP SIZE ADJUSTMENT ' ALGOR ITHM 

GRAD GRADIENT VECTOR GENERATOR 

ITPRT ITERATION OUTPUT GENERATOR 

MINFCN OPTIMI AZTION ALGORITHM 

PINDX PERFORMANCE INDEX GENERATOR 


ELM 

ELMA 

GMADD 

gmprd 

GMSIJR 

GMTRA 

GTOSYM 

LOC 


MCPY 


MM A DO 


M PRO 


MTRA 


SYMTOG 


MATHEMATICAL OPERATIONS PACKAGE 

ALGORITHM TO RETURN VALUE OF AN ELEMENT IN A MATRIX 

ROUTINE TO WRITE INTO AND READ FROM MEMORY AN ELEMENT IN A MATRIX 

ROUTINF TO PERFORM MATRIX ADDITION 

ROUTINE TO PERFORM MATRIX MULTIPLICATION 

ROUTINE TO PERFORM MATRIX SUBTRACTION 

ROUTINE TO 1RANSPOSE A MATRIX, PRESERVING THE ORIGINAL MATRIX 
ROUTINE TO CONVERT A SOUARE N BY N SYMMETRIC MATRIX INTO 
SUPPRESSED SYMMETRIC STORAGE MODE 

ALGORITHM TO COMPUTE A VFCTOR SUBSCRIPT FOR AM ELEMENT IN 
A MATRIX OF SPECIFIED STORAGE MOOE 

ROUTINE TO COPY ONE MATRIX INTO ANOTHER: THE STORAGE MODES 
OE BOTH MATRICES MUST PE IDENTICAL 

ALGORITHM TO FORM THE COMBINATION, C= AL PH A* A+ BF T A* B WHERE 
A,R AND C ARE VECTORS 

ROUTINE TO PERFORM MATRIX MULTIPLICATION; THF TWO MATRICES 
MAY HAVE DIFFERING STORAGE MOOES 

ROUTINE TO TRANSPOSE A MATRIX, PRESERVING THE ORIGINAL MATRIX; 
BOTH MATRICES MUST HAVE THE SAME STORAGE MODE 

ROUTINE TO CONVERT A SYMMETRIC MATRIX (IN SUPPRESSED SYMMETRIC 
STORAGE) INTO A SOUARE N BY N SYMMETRIC MATRIX 


INPUT OUTPUT OPERATIONS PACKAGE 


I MXRNP 

INPUT-OUTPUT 

ROUTINE 

FOR 

IRANDP 

INPUT-OUTPUT 

ROUTINE 

FOR 

MXRNP 

INPUT-OUTPUT 

ROUTINE 

FOR 

NAMRNP 

INPUT-OUTPUT 

ROUTINE 

FOR 

R ANDP 

INPUT-OUTPUT 

ROUTINE 

FOR 

RANDPD 

INPUT-OUTPUT 

ROUTINE 

FOR 


INTEGER-VALUED MATRICES 
INTFGER DATA 
REAL-VALUED MATRICES 
CHARACTER-NAME MATRICES 
CHARACTER DATA 
RFAL DATA 
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35 
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39 

39 

40 

40 
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41 
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43 

44 

44 

45 
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VARIABLE LIST 


EXPERIMENTAL ACTIVE MIRROR PACKAGE 

AIM INVERSE OF AM, TEMPORARY STORAGE 

AM DEFORMATION-FORCE MATRIX, TEMPORARY STORAGE 

ARM REDUCED VERSION OF AM 

ARRM REDUCED VERSION OF ARM 

ASCALV ACTUATOR COMMAND SCALE FACTORS 

DA TEMPORARY STORAGE FLOATING POINT 

DTITCS ACTUATOR SET TIME FOR ITERATIVE CONTROL SYSTEM 

DUMV TEMPORARY STORAGE VECTOR 

FSCALV FIGURE SENSOR MEASUREMENT SCALE FACTORS 

GAINM CONTROL SYSTEM GAIN MATRIX 

GAINV CONTROL SYSTEM GAIN VECTOR 

I INDEX 

IA TEMPORARY STORAGE INTEGER 

ICHNG SWITCH FOR MODIFYING VARIABLE VALUE 
J INDEX 

K INDEX 

L INDEX 

L ACT V ACTUATOR LOCATION DESIGNATOR 
MODE TYPEWRITER CONTROL MODE 
MODES DESIRED VALUE OF MODE 

MODOP DESIRED CONTROL SYSTEM OPERATING MODE 
MODOPS DESIRED VALUE OF MODOP 
N NUMBER OF FIGURE MEASUREMENTS 

NR NUMBER OF ACTUATORS 

NRA DIMENSION OF GAINM 

N SN SWT NUMBER OF SENSE SWITCH USED TO INTERROGATE TYPEWRITER 

TREAL REAL TIME VARIABLF 

UF A V MEASURED ACTUATOR OUTPUTS 

UFV ACTUATOR INPUTS 

XFRV REDUCED XFV VECTOR 

XFSV X COORDINATES OF MEASUREMENT POINTS 
XFV MEASURED FIGURE ERRORS 
XV MAIN DATA BASE 

YFSV Y COORDINATES OF MEASUREMENT POINTS 


PARAMETER OPTIMIZATION PACKAGE 

EPS STEP SIZE 

GRADV GRADIENT VECTOR 

NHC NUMBER OF HALVINGS COMPLETED 

NHM MAXIMUM NUMBER OF STEP SIZE HALVINGS 

NIC NUMBER OF ITERATIONS COMPLETED 

NIM MAXIMUM NUMBER OF SUCCESSFUL ITERATIONS 

NITPRT PRINT OUTPUT EVERY NITPRT ITERATIONS 

NPAR NUMBER OF ADJUSTABLE PARAMETERS 

PARV ADJUSTABLE PARAMETER VECTOR 

PINDEX PERFORMANCE INDEX VALUE 

PISTOR PERFORMANCE INDEX AT END OF LAST SUCCESSFUL ITERATION 
SLGV STORED VALUES OF GRADIENT VECTOR LENGTH 
SPIV STORED VALUES OF PERFORMANCE INDEX 
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EAM PACKAGE 


NAME MAIN 

PURPOSE TO INITIALIZE AND START ACTIVE MIRROR, TO ENGAGE REALTIME OPERATION 

MONITOR, AND TO GENERATE ACTUATOR COMMANDS 

SUBROUTINE REQUIRED INIT 

IRANDP 

MFCS 

MODCHK 

RANDP 

REALT 

TYPCON 

INPUT DATA 
NSNSWT 
( 110 ) 


NAME ACT CAL 

PURPOSE TO TEST FIGURE ACTUATORS 

USAGE CALL ACTCAH NENTRY ) 

PARAMETERS NENTRY= 1 READ IN INITIALIZATION DATA 

NENTR Y=2 ACTUATOR CALIBRATION 

SUBROUTINES REQUIRED ACTCMD 

INIT 

IRANDP 

MXRNP 

RANDPD 

INPUT DATA 
NMEAS 
( 110 ) 

DACT 

( E10.0) 
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NAME ACTCMD 

PURPOSE TO SCALE AND TRANSFER ACTUATOR COMMANDS AND ACTUATOR OUTPUT 

MEASUREMENTS. 

USAGE CALL ACTCMOI NENTR Y ♦ 1 , UF V ,UF AV ♦ ASC AL V ) 

PARAMETERS NENTRY= 1 INITIALIZATION 

NENTRY=2 SCALE AND TRANSFER ACTUATOR COMMANDS 

I DENOTES THE I TH ELEMENT OF UFV AND ASCALV TO BE USED 

UFV ACTUATOR INPUT VECTOR 

UF A V MEASURED ACTUATOR OUTPUT VECTOR 

ASCALV ACTUATOR COMMAND SCALE FACTOR VECTOR 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


CHNG 

TO SET X EQUAL TO V IF I CHNG IS EQUAL TO 1 
CALL CHNGI I CHNG» V » X ) 

I CHNG IF I CHNG =1, SET X= V 

V REAL VARIABLE 

X REAL VARIABLE 
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NAME 


F IGSEN 


PURPOSE TO MEASURE THE FIGURE ERROR AT A DISCRETE POINT ON THE REFLECTING 

SURFACE OF THE MIRROR 

USAGE CALL F I GSEN ( NENTR Y , I , XF V , XF SV , V F SV » F SC ALV ) 

PARAMETERS NENTRY= 1 INITIALIZATION 

NENTRY=2 INTERROGATES FIGURE SENSOR AND SCALES MEASUREMENTS 

I DENOTES THE I TH ELEMENT OF XFSV AND YFSV TO BE USED 

XF V THE FIGURE ERROR VECTOR 

XFSV THE X COORDINATE VECTOR 

YFSV THE Y COORDINATE VECTOR 

FSCALV THE FIGURE ERROR SCALING VECTOR 

SUBROUTINES REQUIRED IRANDP 

RANDPO 

INPUT DATA 
NMEAS 
< 110 ) 

PSCALE 
( E10.0) 


NAME INIT 

PURPOSE TO INITIALIZE THE EXPERIMENTAL ACTIVE MIRROR 

USAGE CALL INITtN ENTRY) 

PARAMETERS NENTRY= 1 TO READ IN INITIALIZATION DATA 

NENTRY=2 INITIALIZE MIRROR FIGURE CONTROL SYSTEM 

SUBROUTINES REQUIRED ACTCMD 

FIGSEN 
IMXRNP 
IRANDP 
MXRNP 
TYPCON 

INPUT DATA 

N NR 

( 2110 ) 

AM 

( RBR 7E10.0) 

FSCALV 
( 7E10.0) 

XFSV 
( 7E10.0) 

YFSV 
( 7E10.0) 

LACTV 
(7110) 

A SC A L V 
( 7E10.0) 
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name 


MFCS 


PURPOSE TO REALIZE MIRROR FIGURE CONTROL SYSTEMS 

USAGE CALL MFCS! NENTRY ) 

PARAMETERS NENTRY= 1 MIRROR FIGURE CONTROL SYSTEM INITIALIZATION 

NENTR Y=2 MEASURE FIGURE ERROR 

NENTRY=3 ACTUATOR OUTPUT CONTROL SYSTEM 

NENTRY=4 ITCS ACTUATOR COMMANDS 

SUBROUTINES REQUIRED ACTCMO 

FIGSEN 

M INFCN 

MPRD 

MTRA 

MXRNP 

RANOPD 

REALT 

REDUAM 

SINV 

INPUT DATA 
GAINV 
( 7E10.0) 

DT ITCS 
I ElO.O) 


NAME MIRCAL 

PURPOSE TO CALIBRATE MIRROR 

USAGE CALL MIRCALI NENTRY) 

PARAMETERS NENTR Y= 1 READ IN INITIALIZATION DATA 

NENTR Y= 2 MIRROR CALIBRATION 

SUBROUTINES REQUIRED ACTCMD 

FIGSEN 

INIT 

IRANDP 

MXRNP 

RANDPD 

INPUT DATA 
NMEAS 
( 110 ) 

DACT 
( ElO.O) 
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NAME 


MODCHK 


PURPOSE 

USAGE 

PARAMETERS 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


NAME 

PURPOSE 

USAGE 

PARAMETERS 

SUBROUTINES 


TO CHECK VALUES OE MODE AND MODOP 

CALL MOOCHKf NENTR Y , MODE , MODED, MODOP , MODOPD , I TEST ) 
NENTRY 

MODE TYPEWRITER CONTROL MODE 

MODED DESIRED VALUE OF MODE 

MODOP DESIRED CONTROL SYSTEM OPERATING MODE 

MODOPD DESIRED VALUE OE MODOP 

ITEST IS SET TO 2 IF MODE DOES NOT = MODED, OR 

IE MODOP DOES NOT = MODOPD 


REALT 

TO INTERROGATE REAL TIME CLOCK 
CALL REALT (TREAL) 

TREAL SUBROUTINE SETS TREAT TO THE REAL TIME 


REDUAM 

TO GENERATE AR AND ARR FROM A 
CALL REDIJAMI NENTRY ) 

NENTR Y= 1 GENERATE AR BY REMOVING COLUMNS FROM A 

NENTR Y=2 GENERATE ARR FROM AR BY REMOVING ROWS FROM AR 

REQUIRED LOC 

MCPY 
MXRNP 
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NAME 


typcon 


PURPOSE 

USAGE 

PARAMETERS 

SUBROUTINES 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


SUBROUTINE 

INPUT DATA 
DPARLV 
( 7E1G.O) 
PARMIN 
C7E10 . 0) 


TYPREWRITER CONTROL OF THE EXPERIMENTAL ACTIVE MIRROR 
CALL TYPCONI NENTRY ) 

NENTR Y= 1 READ IN INITIALIZATION DATA 

NENTRY=? TYPEWRITER CONTROL 

REQUIRED ACTCAL 

CHNG 
ELMA 
IRANDP 
MIRCAL 
REALT 


PARAMETER OPTIMIZATION PACKAGE 


EPCHNG 

TO LIMIT THE ABSOLUTE VALUE OF THE CHANGE IN THE COMPONENTS OF THE 
PARAMETER VECTOR. 

CALL EPCHNGI NENTR Y , STEP , STPDEC , DJDPV ) 

NENTR Y= 1 INITIALIZE DATA 

NENTRY=2 SELECT EPS BASED ON THE ABSOLUTE VALUE OF CHANGE IN 

THE PARAMETER VECTOR. 

STEP STEP SIZE 

STPDEC STEP DECREMENT VALUE 

DJDPV GRADIENT VECTOR 

REQUIRED MXRNP 
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GRAD 


NAME 


PURPOSE 

USAGE 

PARAMETERS 

SUBROUTINES 


TO CALCULATE GRADIENT VECTOR 

CALL GRAD(NENTRY) 

NENTR Y= 1 INITIALIZE DATA 

NENTRY=2 CALCULATE GRADIENT 

REQUIRED ANGRAD 

IRANDP 
PINDX 
RANDPD 


INPUT DATA 
DP AR 
( ElO.O) 
NGRAD 
( 110 ) 


NAME 

ITPRT 


PURPOSE 

TO PRINT 

OPTIMIZATION STATUS 

USAGE 

CALL I TPRT ( NENTR Y ) 

PARAMETERS 

NENTRY= 1 

ITERATION OUTPUT 


NENTRY=2 

RETURN 


NENTR Y=3 

RETURN 


NENTR Y=4 

RETURN 


NENTRY= 5 

ITPRT INITIALIZATION 

SUBROUTINES 

REQUIRED 

DESX 


IRANDP 

MXRNP 

PINOX 


INPUT DATA 
N I TPRT 
( 110 ) 
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NAME 


MINFCN 


PURPOSE TO MINIMIZE A PERFORMANCE INDEX 

OSAGE CALL MI NFCNI NENTRY ) 


PARAMETERS NENTRY= 1 INITIALIZE DATA 

NENTR Y=2 INITIALIZE PERFORMANCE INDEX 

NENTR Y=3 PERFORM MINIMIZATION 


SORROOTINES REQUIRED AVGRAD 

CNGRAD 

DAVIDN 

EPCHNG 

GRAD 

GRADMX 

IRANDP 

ITPRT 

MINFA 

MXRNP 

NEWRAF 

PINDX 

POWEL 

RANDP 

RANDPD 


INPUT DATA 

EPS EPSINC EPSDEC 

( 3E10.0) 

NIM NHM NOPT NPAR ISDEC 

(6110) 

NPAR 

( 110 ) 

PAR V 
( 7E10.0) 


NAME 


PINDX 


PURPOSE TO EVALUATE FIGURE PERFORMANCE INDEX 

USAGE CALL P I NDX ( NENTR Y ) 

PARAMETERS NENTRY= 1 INITIALIZATION 

NENTR Y=2 CALCULATE PINDEX 

SUBROUTINES REQUIRED FIGSEN 


NSDEC 
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MATHEMATICAL OPERATIONS PACKAGE 

NAME ELM 

PURPOSE TO RETURN THE VALUE OF THE L » MTH ELEMENT OF THE MATRIX A 

USAGE X= ELM(A,L,M,N) 

PARAMETERS A A MATRIX WITH N COLUMNS AND AN ARBITRARY NUMBER 

OF ROWS. 

L THE ROW NUMBER OF THE ELEMENT BEING RETURNED 

M THE COLUMN NUMBER OF THE ELEMENT BEING RETURNED 

N THE NUMBER OF COLUMNS IN MATRIX A 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


ELMA 

t 

TO WRITE INTO AND READ FROM MEMORY THE I,JTH ELEMENT OF MATRIX A 
CALL ELMAI NENTRY, A, I , J,V,N) 

NENTR Y= 1 THE VALUE OF V IS STORED INTO THE I,JTH ELEMENT OF A 

NENTRY=2 THE VALUE OF THE I,JTH ELEMENT OF A IS STORED INTO V 

A THE MATRIX WHOSE I,JTH ELEMENT IS USED 

I THE ROW OF THE ELEMENT IN MATRIX A 

J THE COLUMN OF THE ELEMENT IN MATRIX A 

V A SCALAR QUANTITY 

N THE NUMBER OF ROWS IN A 
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NAME 

PURPOSE 

USAGE 

PARAMETERS 


GMADD 

TO PERFORM MATRIX ADDITION, R = A + B 
CALL GMADDI A,B,R,N,M) 

A AN N BY M MATRIX 

B AN N BY M MATRIX 

R AN N BY M MATRIX 

AND R 
B, AND R 


N 

M 


THE NUMBER OF ROWS IN A , B , 
THE NUMBER OF COLUMNS IN A, 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


GMPRD 

TO FORM THE PRODUCT, R = A*B 
CALL GMPRDI A,B,R,N,M,L ) 

A AN N BY M MATRIX 

B AN M BY L MATRIX 

R AN N BY L MATRIX 

N THE NUMBER OF ROWS IN A AND IN R 

M THE NUMBER OF COLUMNS IN A ; ALSO THE NUMBER OF ROWS 

IN B. 

L THE NUMBER OF COLUMNS IN B; ALSO, THE NUMBER OF 

COLUMNS IN R. 
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NAME 


PURPOSE 

USAGE 

PARAMETERS 


GMSUB 

TO PERFORM MATRIX SUBTRACTION, R= A-B 
CALL GMSUBI A,B,R,N,M) 

A AN N BY M MATRIX 

B AN N BY M MATRIX 

R AN N BY M MATRIX 

N THE NUMBER OF ROWS IN A , B , AND R 

M THE NUMBER OF COLUMNS IN A , B , AND R 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


GMTRA 

TO TRANSPOSE MATRIX A INTO MATRIX R 
CALL GMTRAI A,R,N,M) 

A AN N BY M MATRIX 

R AN M BY N MATRIX 

N THE NUMBER OF ROWS IN A: ALSO, THE, NUMBER OF COLUMNS 

IN R. 

M THE NUMBER OF COLUMNS IN A; ALSO, THE NUMBER OF ROWS 

IN R. 
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NAME 


GTOSYM 


PURPOSE TO CONVERT A SQUARE SYMMETRIC NX RY NX MATRIX INTO A SYMMETRIC 

MATRIX IN SUPPRESSED SYMMETRIC STORAGE. THE MATRIX PRODUCED 
CONSISTS OE THE UPPER TRIANGLE (INCLUDING THE DIAGONAL ELEMENTS) 
OF THE NX BY NX MATRIX. THE LENGTH OF THE VECTOR CREATED IS 



NX*< NX + 1 

) /2 . 





USAGE 

CALL GTO 

SYM( X,XS, 

NX ) 




PARAMETERS 

X 

THE 

NX BY 

NX 

SYMMETRIC 

MATRIX (INPUT) 


XS 

THE 

UPPER 

TRIANGLE FORM 

PRODUCED (OUTPUT) 


NX 

THE 

ORDER 

OF 

THE INPUT 

MATRIX 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


LOC 

TO COMPUTE 
SPECIFIED 

CALL LOCI I 

I 

J 

IR 

N 

M 

MS 


A VECTOR SUBSCRIPT FOR AN ELEMENT IN A MATRIX OF 
STORAGE MODE. 

, Jt IR*N»M,MS) 

ROW NUMBER OF ELEMENT 

COLUMN NUMBER OF ELEMENT 

RESULTANT VECTOR SUBSCRIPT (DETERMINED BY LOC) 
NUMBER OF ROWS IN MATRIX 
NUMBER OF COLUMNS IN MATRIX 

A SINGLE DIGIT INDICATING THE STORAGE MODE OF THE 
MATRIX. 

0 GENERAL 

1 SYMMETRIC 

2 DIAGONAL 
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NAME 


MCPY 


PURPOSE TO COPY MATRIX A INTO MATRIX R 

USAGE CALL MCP Y< A , R ,N , M , MS ) 

PARAMETERS A AN N BY M MATRIX 

R AN N BY M MATRIX 

N THE NUMBER OF ROWS IN A AND IN R 

M THE NUMBER OF COLUMNS IN A AND IN R 

MS A SINGLE DIGIT INDICATING STORAGE MODE OF ROTH A 

AND R 

0 GENERAL 

1 SYMMETRIC 

2 01 AGONAL 

SUBROUTINES REQUIRED LOC 


NAME 

MM ADD 



PURPOSE 

TO FORM THE 

COMBINATION, C = ALPHA*A + BETA<=R 

USAGE 

CALL MMADOI N 

, ALPHA, 

BETA, R,C ) 

PARAMETERS 

N 

LENGTH 

OF VECTORS A,B, AND C 


ALPHA 

SCALAR 

QUANTITY 


A 

VECTOR 

OF LENGTH N 


BETA 

SCALAR 

QUANTITY 


R 

VECTOR 

OF LENGTH N 


C 

VECTOR 

OF LENGTH N 
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NAME 

PURPOSE 

USAGE 

PARAMETERS 


SUBROUTINES 

NAME 

PURPOSE 

USAGE 

PARAMETERS 


SUBROUTINES 


MPRD 


TO FORM THE PRODUCT, R=A*B 


CALL MPRO( A,B,R,N,M,MSA,MSB»L) 


A 

B 

R 

N 

M 


MSA 


MSB 


L 


AN N BY M MATRIX 
AN M BY L MATRIX 
AN N BY L MATRIX 

THE NUMBER OF ROWS IN A AND IN R 

THE NUMBER OF COLUMNS IN A? ALSO, THE NUMBER OF ROWS 
IN B. 

A SINGLE DIGIT INDICATING STORAGE MODE FOR MATRIX A 

0 GENERAL 

1 SYMMETRIC 

2 DIAGONAL 

A SINGLE DIGIT INDICATING STORAGE MODE FOR MATRIX R 

0 GENERAL 

1 SYMMETRIC 

2 DIAGONAL 

THE NUMBER OF COLUMNS IN B; ALSO, THE NUMBER OF 
COLUMNS IN R. 


REQUIRED LOC 


MTRA 


TO TRANSPOSE MATRIX A TO FORM MATRIX R 
CALL MTRA( A,R,N,M,MS) 

A AN N BY M MATRIX 

R AN M BY N MATRIX 

N THE NUMBER OF ROWS IN A; ALSO, THE NUMBER OF COLUMNS 

IN R. 

M THE NUMBER OF COLUMNS IN A; ALSO, THE NUMBER OF ROWS 

IN R. 

MS A SINGLE DIGIT INDICATING THE STORAGE MODE OF BOTH A 

AND R. 

0 GENERAL 

1 SYMMETRIC 

2 DIAGONAL 

REQUIRED MCPY 
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NAME 


SYMTOG 


PURPOSE TO CONVERT A SYMMETRIC MATRIX (IN SUPPRESSED SYMMETRIC STORAGE), 

WHOSE LENGTH IS NX*( NX + 1 ) /2 , INTO A SQUARE SYMMETRIC NX BY NX 
MATRIX WHOSE LENGTH IS NX*NX. 

USAGE CALL S YMTOGt X S , X , NX ) 

PARAMETERS XS THE SYMMETRIC MATRIX VECTOR (INPUT) 

X THE EXPANDED GENERAL MATRIX VECTOR (OUTPUT) 

NX THE ORDER OF THE INPUT MATRIX 


INPUT OUTPUT OPERATIONS PACKAGE 

IMXRNP 

TO READ AND PRINT I NT EGER-V ALIJED MATRICES 
CALL IMXRNP( MA,NA,NB,NENTRY) 

MA 1-DIMENSION INTEGER VECTOR WHOSE DIMENSION MUST BE AT LEAST 
NA*NB. MATRIX IS STORED COLUMN-WISE IN THIS VECTOR, 

NA NUMBER OF ROWS IN MATRIX 

NB NUMBER OF COLUMNS IN MATRIX 

NENTR Y= 1 READ IN AND PRINT OUT MATRIX 

NENTRY=2 READ IN MATRIX 

NENTRY=3 PRINT OUT MATRIX 

NENTR Y=4 READ IN HEADING CARD, READ IN MATRIX AND 

PRINT OUT MATRIX. 

NENTRY=5 PUNCH OUT MATRIX 

REQUIRED RANDP 

INTEGER MATRIX (ROW-WISE) 

PRINTED INTEGER MATRIX (ROW-WISE) 


NAME 

PURPOSE 

USAGE 

PARAMETERS 


SUBROUTINES 
INPUT DATA 
OUTPUT 
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NAME 


IRANDP 


PURPOSE 

USAGE 

PARAMETERS 


SURROUTINES 
INPUT DATA 

OUTPUT MAY 

NAME 

PURPOSE 

USAGE 

PARAMETERS 


SURROUTINES 
INPUT OATA 
OUTPUT 


TO REAO AND PRINT INTEGER OATA 


CALL IRANDPI NO,IA,IR,IC,ID,IE,IF,IG,NENTRY) 


ND 

THE NUMBER OF INTEGER VALUES TO BE READ IN 
PRINTED OUT. 

AND 

I A- 1 G 

THE INTEGER VARIABLES TO WHICH THE INTEGER 
WILL BE ASSIGNED. 

values 

NENTRY= 1 

READ IN 7 INTEGER VALUES (7110) FOR THE INTEGER 
VARIABLES IA-IG AND PRINT THE FIRST ND VARIABLES. 

NENTRY=2 

PRINT THE FIRST ND INTEGER VARIABLES 


NENTRY=3 

RETURN 


NENTR Y=4 

READ AND PRINT HEADING CARD, READ 7 INTEGER 
THE FIRST ND VARIABLES. 

(7110) FOR THE INTEGER VARIABLES IA-IG AND 

VALUES 
PRINT OUT 

REQUIRED 

RANDP 


MAY BE 

HEADING CARD 
INTEGER VALUES CARD 


BE 

PRINTED HEADING 
PRINTED INTEGER VALUES 



MXRNP 

TO READ AND PRINT REAL-VALUED MATRICES 
CALL MXRNPI VA , N A , NB , NENTR Y > 

VA 1-DIMENSION REAL VECTOR WHOSE DIMENSION MUST BE AT LEAST 
NA*NB. MATRIX IS STORED COLUMN-WISE IN THIS VECTOR. 

NA NUMBER OF ROWS IN MATRIX 

NB NUMBER OF COLUMNS IN MATRIX 

NENTRY= 1 READ IN AND PRINT OUT MATRIX 

REAO IN MATRIX 

PRINT OUT MATRIX 


NENTR Y=2 
NENTRY=3 
NENTR Y=4 

NENTR Y= 5 
REQUIRED 


READ IN HEADING CARD, READ IN MATRIX AND 
PRINT OUT MATRIX. 

PUNCH OUT MATRIX 

RANDP 

REAL MATRIX ( ROW-WISE ) 

PRINTED REAL MATRIX (ROW-WISE) 
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NAME 

NAMRNP 


PURPOSE 

TO READ, PRINT AND STORE MATRIX M WHICH CONTAINS 4-CHARACTER NAMES 

USAGE 

CALL NAMRNPI M,NA,NR,NENTRY ) 

PARAMETERS 

M 

AN NA BY NB MATRIX 


NA 

THE NUMBER OF ROWS IN M 


NB 

THE NUMBER OF COLUMNS IN M 


NENTRY= 1 

READ IN AND PRINT M AND STORED M INTO A SINGLE- 
DIMENSIONED VECTOR. 


NENTRY=2 

READ IN MATRIX M AND STORE M INTO A SINGLE- 
DIMENSIONED VECTOR. 


NENTRY=3 

PRINT OUT MATRIX M 


N ENTR Y = 4 

READ IN A HEADING CARD, READ IN MATRIX M AND 

STORE M INTO A S I NGLE-D 1 MENS I ONED VECTOR, AND PRINT M 

SUBROUTINES 

REQUIRED 

RANDP 

INPUT DATA 

AN NA BY 

NB MATRIX OF 4-CHARACTER NAMES 

OUTPUT 

PRINT THE 

NA BY NB MATRIX 


NAME 

RANDP 


PURPOSE 

TO READ AND 

PRINT HEADING CARDS 

USAGE 

CALL RANDPI NENTRY ) 

PARAMETERS 

NENTRY= I 

READS CARD IN BAB FORMAT AND PRINTS CONTENTS (AT TOP OF 
NEXT PAGE) IN BAB FORMAT. 


NENTRY=2 

READS CARD IN 8A8 FORMAT AND PRINTS CONTENTS IN 
BAB FORMAT. 


NENTRY=3 

READS CARD IN 7(2X,AR) FORMAT AND PRINTS CONTENTS IN 
7 ( 7 X , A B ) FORMAT AT TOP OF NEXT PAGE. 


NENTR Y=4 

READS CARD IN 7(2X,A8) FORMAT AND PRINTS CONTENTS IN 
7(7X,A8) FORMAT. 

INPUT DATA 

HEADING CARD 


OUTPUT 

PRINTED HEADING 
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NAME 


RANDPD 


PURPOSE TO READ AND PRINT FLOATING POINT DATA 

USAGE CALL R ANDPDI ND , DA , DB , DC , DE , DF , DG,NENTR Y ) 

PARAMETERS ND THE NUMRER OF REAL FLOATING PT. VALUES TO RE READ IN 

AND/OR PRINTED OUT. 

DA-DG THE REAL VAR1ARLES TO WHICH THE REAL VALUES WILL BE 

ASSIGNED. 

NENTR Y= 1 READ IN 7 REAL VALUES (7EIO.O) FOR REAL VAR I ARLES 
DA-DG AND PRINT THE FIRST ND VARIABLES. 

NENTRY=2 PRINT OUT THE FIRST ND VARIABLES. 

NENTR Y=3 READ IN REAL VALUES (7E10.0) FOR REAL VARIABLES DA-DG 

NENTRY=4 READ AND PRINT HEADING CARD, READ IN 7 REAL VALUES 

(7E10.0) FOR THE REAL VARIABLES DA-DG, AND PRINT OUT 
THE FIRST ND VARIABLES. 

SUBROUTINES REQUIRED RANDP 

INPUT DATA MAY BE HEADING CARD 

FLOATING PT. VALUES CARD 

OUTPUT MAY BE PRINTED HEADING 

PRINTED FLOATING PT. VALUES 
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CHAPTER 5 

TYPICAL DATA CONFIGURATIONS 


5.0 Introduction 


Data for the EAM is introduced in two different fashions: card input 
and typewriter input. Typewriter input data is detailed in section 3. 2. 
Card input data comprises the bulk of the information utilized by the 
EAM software. The following sections discuss the card data deck 
structure 


5.1 Card Input Data for the SLCS and LOCS 


If Simplified Linear ( SLCS ) or Linear Optimal ( LOCS ) Control is 
desired, the following data deck is read by the card reader: 


NSNSWT 

( 2110 ) 

N NR 

( 2110 ) 

AM 

(RBR 7E10. 0) 

FSCALV 
(7E10. 0) 

XFSV 
(7E10. 0) 

YFSV 
(7E10. 0) 

LACTV 
( 7E10. 0) 

ASCALV 
( 7E10. 0) 

NMEAS 

(110) 

PSCALE 
(E10. 0) 

NMEAS 

(110) 

DACT 
(E10. 0) 

NMEAS 

( 110 ) 

DACT 
(E10. 0) 

GAINV 
(7E10. 0) 

DTITCS 4 

(E10. 0) 



The data deck includes variable heading cards, to assist identification, 
as well as the numerical data in the indicated format. The heading 
card is read and printed by the computer over the corresponding data 
in the output print. The notation (RBR 7E10. 0) under the heading AM 
indicates that the variable AM (the force deformation matrix of the 
mirror) is read in row by row in 7E10. 0 (unassigned decimal) format. 
Each row must be started on a new card. 

5.2 Card Input Data for ITCS 

If the iterative figure control system (ITCS) is to be realized, the soft- 
ware will read in the data deck described in section 5.1, and the following 
set of heading and numerical data cards: 


EPS EPSINC EPSDEC 

(3E10. 0) 

NIM NHM NOPT 

(6110) 

NPAR 

( 110 ) 

PARV 
(7E10. 0) 

DPAR 
(E10. 0) 

NGRAD 

( 110 ) 

DPARLV 
(7E10. 0) 

PARMIN 
(E10. 0) 

NITPRT 

( 110 ) 


NPAR IS DEC NSDEC 
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CHAPTER 6 

SPECIFICATION OF A SMALL HYBRID COMPUTER FOR THE EAM 


6„ 0 Introduction 


A key part of the study of advanced space telescope technology currently 
in progress at NASA-MSFC and MIT-DL will be an experimental program 
which will lead to the evaluation and development of techniques for ac- 
tively controlling the optical surface shape or figure of primary telescope 
mirrors,, A necessary feature of the experimental program is a control 
system which will feature an element (or set of elements) which are 
easily and inexpensively modified to reflect control algorithm design changes. 

Adopting a modern approach, it was decided to use a small general 
purpose digital computer as the programmable control system component. 
Such an approach is economical as a result of the significant reduction 
in small computer manufacturing costs, which have occurred over the 
last few years, and the high cost of constructing special purpose analog 
systems. 

The digital computer also provides a valuable tool for real-time processing 
of experimental data and the solution of other control problems associated 
with system alignment. 

The main disadvantage of the digital computer is the character of the 
information it handles. Each piece of information in the computer is 
characterized by a binary number. The number is expressed by one 
dimensional array of N elements each of which can have two values 
which are identified by 0 and 1. A power of 2 is associated with each 
element. Hardware devices in control systems, on the other hand, 
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generally produce data which is of a one dimensional, or analog, character. 
Thus conversion devices must be used which map the one dimensional 
analog signals into N dimensional binary numbers for computation purposes. 
Converters are also required to map N dimensional binary data into one 
dimensional signals for use by hardware elements. Both of these oper- 
ations are easily performed using high-speed digital to analog and analog 
to digital converters. 

After discussions with a number of people at MIT-DL who are using 
the Digital Equipment Corporation PDP-9 (a predecessonof the PDP-15) 
and Dean Hamilton (formerly in charge of the PDP-10 digital computer 
facility at NASA-ERC and currently director of the CARS computer at 
MIT-DL), it was decided to restrict attention to the Digital Equipment 
Corporation PDP-15 Computer. 12 12 

6. 1 A Small Computer for the EAM 

6,1.1 Manufacturer 

Digital Equipment Corporation 
146 Main Street 

Maynard, Massachusetts 01754 

6. 1. 2 Manufacturer' s Designation 
Model PDP-15 
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G. 1. 3 Basic Computer 


The basic computer consists of a central processor and a random access 
core memory. The central processor is capable of executing a set of 
instructions which include fixed point arithmetic operations on data (add, 
subtract, multiply and divide); data transfer operations to and from 
the core memory, mass storage devices, and peripherals such as the 

s'c 

teleprinter and other data displays; and logic operations such as test 
and branch. The basic core memory consists of 4, 096 words. The 

3*C jfc 

memory cycle time is 0. 8 microseconds which is relatively fast. 

6. 1. 4 Memory Expansions 

The random access core memory can be expanded in 4, 096 word incre- 

* * * 

ments to a maximum of 32, 768 words. 


Tape units, disks 

Cathode ray plotter, calcomp plotter, high-speed printer 

Memory "Cycle time" is defined as the time required to trans- 
fer a data word to or from the core memory. 

Further expansion to 131, 072 words is possible with hardware 
modifications. 
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6. 1. 5 Analog to Digital Conversion 


The conversion is performed by one converter which is sequentially 
switched to each input channel. Each channel level is converted to a 
binary number which is transferred to the central processor. The 
analog signal switching is performed by field effect transistor switches. 
The converter features a selectable 6 to 12 bit output word length which 
provides more than enough resolution for handling figure sensor data. 

6. 1. 6 Digital to Analog Conversion 

Digital to analog conversion is normally performed by transferring the 

>;< sjc 

binary number into a buffer register where it is stored. The stored 
binary number is then observed by a device which produces an analog 
signal proportional to the value of the binary number. The resolution 
in the analog signal is determined by the number of bits in the buffer 
register. A length of at least 10 bits is required to yield acceptable 
resolution (1 part in 1,024). 


* If the maximum figure error is 5\ and the minimum sensed error is 
A/200, the ratio (5A)/(A/200) = 1, 000 <2 1U = 1,024. The actuator 
commands are linearly related to the figure errors: thus, the 
same criteria may be applied to the output signals from the computer 
which drives the actuators. 

** An operational amplifier with an appropriate network of input 
resistors. 
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6.1.7 Mass Memory 

DEC manufactures a wide variety of tape transports and disc memories 
for the PDP-15. The small Dectape units are inexpensive and can 
virtually eliminate the necessity of using paper tape or cards for storage 
(a single reel of Dectape can store 150, 000 18 bit words). 

6. 1. 8 Data Communications 

The PDP-15 is designed for use with a variety of teletypes, high-speed 
printers, paper-tape readers and punches, card readers and cathode 
ray plotting devices. 

6.1.9 Special Options 

A particularly useful control processor option is the extended arithmetic 
element which reduces the time required to perform multiplications 
and divisions by a factor of approximately 50. Since computer utili- 
zation in a real-time environment is usually limited by computation speed 
rather than memory capacity, reducing the time required to perform 
multiplications and divisions (the most time-consuming operations) can 
significantly increase the usefulness of the computer. 

6. 2 Component Cost for PDP-15 

6. 2. 1 Basic System 

Central Processor 

4, 096 18 bit words of core memory PDP-15 $15, 600 
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6. 2.2 Memory Expansion for Basic System 

4,096 word core module to increase 
Basic System memory capacity to 

8, 192 words MK15-A 

6. 2. 3 Memory Modules for Further Expansion 
Assembly to further increase memory 

capacity by 8, 192 words MM15-A/MK15-A 

6. 2. 4 Console Teleprinter (Heavy Duty) KSR-35 

6. 2. 5 High-Speed Paper Tape Reader 

and Punch PC 15 

Paper tape input is required to run 
diagnostic routines 

6. 2. 6 Analog to Digital Conversion 

6 to 12 bit selectable analog to 
digital converter and multiplexer 

control AF03B 

9-35 pi sec conversion time 

multiplexes up to 64 channels A121 

4-Channel field effect transistor 
switch: One required for each set 
of four analog to digital channels 

Scaling amplifiers (if required) AF01B 

Convert input signals to an average 


6,000 

14,000 

3,000 

3, 900 

5, 000 
65 

300 
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level and amplitude suitable for 
the analog to digital converter and 
FET switches (one for each channel) 

6. 2. 7 Digital to Analog Conversion 

Multiplexer control for up to 16, 10 bit. 


digital to analog converter channels 

AA05A 

$ 5,500 

Expansion of AA05A to 64 channels 

AA05B 

2, 900 

Digital to analog converter, single 
buffered ± 5. 0 volts or ± 10. 0 volts: 
One required per channel 

A609 

375 

Reference power supply: For up to 
12 digital to analog converters 

A 610 or A611 

400 

6. 2. 8 DEC Tape Drives 



DEC tape control unit for 8 tape units 

TC02D 

5, 400 

Dectape transport 

TU55 

2, 350 

6. 2. 9 Extended Arithmetic Element 

KE1 5 

2, 800 

Decreases multiply and divide time 
by approximately 50 



6. 2. 10 Positive to Negative Buss Converter 

DW15A. 

2, 000 

Required for compatibility by some 
D/A and A/D modules 



6.2.11 Disk Memory 

RP02 

9, 000 
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Stores up to 262, 144 words per 
disk, average access time 16.7 
milliseconds 

Control unit for up to 8 RP02 disks RP15 $6, 000 

6. 3 Prices of Typical System Configurations 

6. 3. 1 Basic System Expanded to 16, 384 Words of Core Memory, 

Teleprinter, Paper-Tape Reader and Punch 


1 (PDP-15/10) 

$15, 600 

2 (MK15-A) 

12,000 

1 (MM15-A) 

8, 000 

1 (KSR-35) 

3,000 

1 (PC 15) 

3,900 

TOTAL 

$42, 500 


6. 3. 2 System 6.3.1 Expanded to Include Three DEC Tape Units 


1 (6. 3. 1) 

$42, 500 

1 (TC02D) 

5, 400 

3 (TU55) 

7,050 

TOTAL 

$54, 950 


Time required to acquire a piece of information stored on the disk: 
Access time is primarily a function of disk angular velocity. 
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6„ 3. 3 Basic Digital to Analog Package (16 Channels) 


1 (AA05A) 

$5, 500 

16 (A 609) 

6, 000 

2 (A 61 0) 

1, 600 

TOTAL 

$13, 100 

Basic Analog to Digital Package (16 Channels) 

1 (AF01B) 

$ 5,000 

4 ( A 1 2 1 ) 

260 

16 (AH03) 

4, 800 

TOTAL 

$10, 060 

System 4. 2 with D/A and A/D Packages 


4. 3 and 4. 4 


1 (3. 2) 

$54, 950 

1 (3. 3) 

13, 100 

1 (3. 4) 

10, 060 

TOTAL 

$78, 110 


6. 4 Summary 

The system defined in Section 3. 5 represents a minimal system for an 
active mirror experimental program. Such a hybrid computer would 
permit configurations using up to 16 actuators and sensor outputs to be 
investigated. More complex systems will probably require an increase 
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in core memory and modifications in the A/D and D/A channeling to 
carry the increased load. 
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CHAPTER 7 

SUMMARY AND CONCLUSIONS 

The software requirements for an experimental active mirror have 
been defined and a complete preliminary software package design has 
evolved. Although the software was specifically designed for com- 
patibility with the XDS Sigma 5/7 system, its FORTRAN II- IV structure 
permits it to be executed by practically any FORTRAN operating system 
with minor modifications. 

The software package realizes the linear optimal, simplified linear 
and iterative control algorithms discussed in reference 1. In addition 
the software provides servo loop closures about the figure actuators 
which help to eliminate the effect of actuator characteristics variations. 

The software package includes routines for actuator calibration and the 
generation of discrete mirror models. An operating package is also 
incorporated which enables the investigator to control the operation of 
the EAM and to operate the system in diagnostic or system parameter 
modification modes. 

It should be emphasized that the software has not been thoroughly debugged. 
Rigorous testing will be an important segment of the follow-on work asso- 
ciated with this project. Testing will include the definition of software 
models of the hardware components so that the software package can 
be completely debugged without the necessity of communicating with 
hardware elements. This will enable most of the software tests to be 
performed using the MIT computer facility. 
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The report concludes with a description of a small hybrid computer 
(DEC PDP 15) for application in an experimental active mirror. While 
the specified computer facility is not directly applicable to the EAM 
as it is currently envisioned it is useful to review the technical consider- 
ations involved in the selection of the PDP 15 and its peripherals. 
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APPENDIX A 
EAM PACKAGE 


A. 0 Introduction 


This appendix contains listings of the following subroutines. 


MAIN PROGRAM 62 
SUBROUTINE ACTCAL ( NENTRY ) 64 
SUBROUTINE ACTCMDI NENTR Y , I , UFV , UE AV , A SC AL V ) 65 
SUBROUTINE CHNGI I CHNG , V , X ) 65 
SUBROUTINE F I G SEN! NENTRY , I , X F V , XF SV , V FX V , F SC AL V ) 66 
SUBROUTINE I N I T ( NENTRY ) 67 
SUBROUTINE MFCS! NENTRY) 68 
SUBROUTINE M I RCAL ( NENTRY ) 71 
SUBROUTINE MODCHKI NENTRY, MODE, MOOED, MODOP , MODOPD , I TEST ) 72 
SUBROUTINE REALTI NENTRY ) 72 
SUBROUTINE REDUAMI NENTRY ) 73 
SUBROUTINE TYPCONI NENTRY ) 74 


61 



c 

MAIN 

EAM10000 

c 


EAM10010 

c 

MAIN PROGRAM FOR EXPERIMENTAL ACTIVE MIRROR 

E AM10020 

c 


EAM10030 


DIMENSION XVI 10000) , AMI 3000) , AIM! 3000) ,XFV( 20 ) , UFV I 20 ) , ASCAL V I 20 ) 

, E AM10040 

1 

FSCALVI 20) ,XFSV( 20) ,YFSV( 20) , XFRVI 20) ,DUMV< 20) ,UFAV< 20) , DUMVA I 20 ) E AMI 0050 

2 

,GAINV( 10) tGAINMI 400) ,LACTV( 20) 

E AMI 0060 


COMMON/BLKMFC/ N , NR , XV , MODE , LACTV , MODOP , NSNSWT * NR A 

EAM10070 


EOIJ I VALENCE I X Vll) , AMI 1 > ) , I XVI 3001 ) » A I MI 1 ) ) , I XV ( 6001 ) , XFV ( 1 ) ) , 

EAM100R0 

1 

I UFVI 1) »XV( 6101) ) ,( ASCAL VI 1 ) , X VI 6 1 2 1 ) ) , I FSCALVI 1 ) ,XV< 6141 ) ) , 

EAM10090 

? 

I XFSVI 1) ,XV< 6161 ) ) ,1 YFSVI 1) ,XV! 6181) >,( XFRVI 1 >»XV< 6201) ) , 

EAM10100 

3 

I UFAVI 1 ) ,XV( 6 221 ) ) ,< DUMVI 1 ) , X VI 6241 ) ) ,( DUMVAI 1) ,XV< 6261 ) ) , 

EAM101 10 

4 

I GAINVI 1 ) tXVI 6281 ) ) ,( GAINMI 1)1 XVI 6291) ) 

EAM10120 

c 


E AMI 0 1 30 

1000 

FORMAT! 1H 1 , / , 4HM A I N ) 

EAM10140 

1010 

FORMAT! 9HMDD0P ERR) 

EAM10150 

c 


EAM10160 

c 

INITIALIZATION 

EAM10170 


PRINT 1000 

EAM10180 


CALL RANDPI 2 ) 

EAM10190 


CALL I RANDPI 1, NSNSWT, IA,IA,IA,IA,IA,IA,4) 

EAM10200 


CALL TYPCONI 1) 

EAM102 10 


CALL INIT(l) 

EAM10220 

2105 

CALL TYPCONI 2) 

EAM102 30 

r. 


EAM10240 

c 

IS CONTROL CONFIGURATION DEFINED 

EAM10250 


IF( MODOP- 5) 21 10,21 10, 2120 

EAM10260 

c 

NO 

EAM10270 

2110 

TYPE 1010 

EAM10280 


GO TO 2105 

EAM10290 

c 

YES 

EAM10300 

2120 

MODOP S=M0DOP 

EAM10310 

C 


EAM10320 

C 

INITIALIZE ACTIVE MIRROR 

E AM10330 


CALL TYPCONI 2) 

E AMI 0340 


MODE S= 1 

EAM10350 


CALL MODCHKI 1 , MODE t MODE S , MODOP , MODOP S , I A ) 

EAM10360 


GO TO I 2130,2105) , I A 

EAM10370 

2130 

CALL INITI2) 

EAM10380 


CALL MFCS! 1 ) 

EAM10390 


IF (MODOP-8) 2133, 2132, 2133 

EAM10391 

2132 

CALL MINFCN Cl) 

EAM10 39 2 


CALL MINFCN (2) 

EAM10393 

c 


EAM10400 

c 

START ACTIVE MIRROR 

EAM10410 

2133 

MODE S=2 

EAM10420 

2135 

CALL TYPCONI 2) 

EAM10430 


CALL MODCHKI 1 , MODE , MODES , MODOP , MODOP S , I A ) 

EAM10440 


GO TO I 2140,2135) ,IA 

EAM10450 

2140 

CALL REALTITREAL) 

EAM10460 


TSTOR E = TRE AL 

EAM10470 


GO TO 2150 

EAM10480 

C 


EAM10490 

C 

RFALTIME OPERATION MONITOR 

E AM10500 

2141 

XVI 61511=0.0 

EAM10510 

2142 

CALL REALTITREAL) 

E AMI 0520 


IF! TREAL-TST0RE-DT)2142, 2143, 2143 

EAM10530 

2143 

TST0RE= TREAL 

EAM10540 

r. 


EAM10550 

C 


E AMI 05 60 
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o n 


C GENERATE ACTUATOR COMMANDS EAM10570 

2150 CALL MFCS( 2 ) EAM10580 

EAM10590 

INTERROGATE TYPEWRITER IF SENSE SWITCH NSNSWT IS SET EAM10600 

2145 I F ( SENSE SWITCH NSNSWT) 2205,2147 EAM10610 

2205 CALL REALTI TREAL ) EAM10620 

0A=TREAL-T STORE EAM10630 

2146 CALL TYPCONI 2 ) EAM10640 

CALL MODCHKI 1, MODE, MOOES, MODOP,MODOPS,IA) EAM10650 

GO TO ( 2147,2146) I A EAM10660 

2147 CALL REALT ( TREAL ) EAM10670 

TSTORE= TREAL -DA EAM10680 

GO TO 2141 EAM10690 

c EAM10700 

END EAM10710 
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ooo o ouu ou uu 


SUBROUTINE ACTCAL I NENTRY ) 

SUBROUTINE TO TEST EIGURE ACTUATORS 

DIMENSION XV( 10000) , AM! 3000) , AIM! 3000) ,XFV( 20 ) , UFV ( 20 ) , ASCALV I 20 ) 

1 FSCALVI 20) ,XFSV( 20) ,YFSV< 20) ,XFRV< 20) , DUMVI 20) ,UFAV( 20) ,DUMVA< 20 

2 , GAINVI 10) ,GAINM< 400) ,LACTV< 20) 

COMMON /BLKMFC/ N , NR , XV , MODE , L AC TV , M000P , NSNSWT , NR A 
EQUIVALENCE ( XVI 1 ) , AM( 1 ) ) , ( XVI 3001 > , AIM! 1 > ) ,< X V I 600 1 ) , XF V I 1 ) ) , 

1 ( UFVI 1) ,XV( 6101 ) ) , l ASCALVI 1 ) , X VI 61 2 1 > ) , I FSCALVI 1 ) ,XV< 6141 ) ) , 

2 I XFSVI 1) ,XV( 6161 ) ) ,< YFSVI 1 ) ,XV! 6181 ) ) , I XFRVI 1 ), XV I 6201 ) ) , 

3 I UF AVI 1 ) , XVI 6 221 ) ) ,( DUMVI 1 ) , X V I 6241 ) ) ,( DUMVAI 1 ) , XV I 6261 ) ) , 

4 ( GAINVI 1 ) f XVI 6281 ) ) , I GA I NMI 1 ) I X VI 6291 ) ) 

1000 FORMAT! 6HACTCAL ) 

1002 FORMAT! 1 HI /3X , 12HACT0UT/ ACT I N ) 

1003 FORMAT! 7X,8HACAL END) 

GO TO I 1,2), NENTRY 


INPUT DATA 
PRINT 1000 
CALL IRANDPI 
CALL RANDPDI 
0B=DACT*2.0 
RETURN 


1 ,NMEAS, IA,IA,IA,IA,IA,IA,4) 
1,DACT,DA,DA,0A,0A,DA,0A,4) 


2201 

2202 


ACTUATOR CALIBRATION 
CALL INITI2) 

DO 2202 1 = 1, NR 
DO 2201 J= 1 , NME AS 
UFVI I ) =-DACT 

CALL ACTCMDI 2 , I , UF V , UF AV , A SC ALV ) 
OA = UF AVI I ) 

UFVI I ) = DACT 

CALL ACTCMDI 2 , I ,UFV ,UFAV , ASCALV ) 
OUMVI I ) = ( UF AVI I ) -DA ) /DB+DUMVI I ) 
DUMVI I ) = DUMVI I ) /NMEAS 

PRINT OUT ACTUATOR SCALE VECTOR 
PRINT 1002 

CALL MXRNPI DU MV, 1 »NR,3) 

TYPE 1003 
RETURN 

END 


E AM10720 
EAM10730 
EAM10740 
EAM10750 
, EAM10760 
) EAM10770 
EAM10780 
EAM10790 
EAM10800 
EAM108 10 
EAM10820 
E AM10830 
EAM10840 
EAM10850 
EAM10860 
EAM10870 
EAM10880 
EAM10890 
EAM10900 
EAM10910 
EAM10920 
EAM10930 
EAM10940 
EAM10950 
E AM 1 09 60 
EAM10970 
EAM10980 
EAM10990 
EAM11000 
EAM11010 
E AMI 1020 
EAM11030 
EAM1 1040 
EAM11050 
E AMI 1060 
EAM11070 
EAM11080 
EAM11090 
E AMI 1 100 
EAM11110 
EAM1 1120 
EAM1 1130 
E AMI 1 140 
EAM11150 
E AMI 1 1 60 
EAM11170 
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SUBROUTINE ACTCMD< NENTRY, I , UFV , UF AV , ASC AL V ) 

E AM 1 1 1 80 

c 


EAM11190 

c 

SUBROUTINE TO SCALE AND TRANSFER ACTUATOR COMMANDS 

AND ACTUATOR 0UEAM11200 

c 

MEASUREMENTS. 

EAM11210 

c 


E AMI 1 220 


DIMENSION UF V( 1 ) ,UFAV( 1 ) ,ASCALV(1) 

E AMI 12 30 

c 


EAM11240 


GO TO( 1,2) ,NENTRY 

EAM11250 

c 


E AMI 1 260 

c 

INITIALIZATION 

E AMI 12 70 

1 

RETURN 

EAM11280 

c 


EAM11290 

c 

SCALE AND TRANSFER ACTUATOR COMMANDS 

EAM11300 

2 

DA = UF V ( I ) * A S C A L V ( I ) 

EAM11310 

c 


EAM11320 

c 


EAM11330 

c 

INSERT DTOA SOFTWARE HERE 

EAM11340 

c 


EAM11350 

c 


EAM11360 

c 

TRANSFER ACTUATOR OUTPUT MEASUREMENTS 

EAM11370 

c 


EAM113R0 

c 


E AM 1 1390 

c 

INSERT ATOD SOFTWARE HERE 

E AM 1 1 400 

c 


EAM11410 

c 


EAM11420 


UEAVl I ) = DA 

EAM11430 


RETURN 

EAM1 1440 

c 


EAM11450 


END 

EAM11460 



SUBROUTINE CHNG( I CHNG, V ,X ) 

E AMI 1470 

C 


E AM 1 1 480 

C 

SUBROUTINE SETS X EQUAL TO V IF I CHNG IS EQUAL TO 1 

E AMI 1490 

c 


EAM11500 


GO TO (1,2) , I CHNG 

E AMI 1510 

c 


EAM11520 

1 

x=v 

EAM11530 

2 

RETURN 

EAM1 1540 

C 


EAM11550 


END 

EAM11560 
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nooooo ooo no o o non 


SUBROUTINE F I G SEN ( NENTRY,I,XFV,XFSV,YFSV,FSCALV) 


EAM11570 


C EAM1 

SUBROUTINE TO MEASURE THE FIGURE ERROR XFVII) AT A DISCRETE POINT E AMI 
COORDINATES XFSVI I > , YFSVII) ON THE REFLECTING SURFACE OF THE MI RROEAM1 

EAM1 

DIMENSION XFV( 1 ) ,XFSV( 1 ) , YFSV< 1 ) EAM1 

E AMI 

1000 FORMAT! 6HFIGSEN) EAM1 

EAM1 

GO T0( 1,2) , NENTRY E AMI 

E AMI 

INITIALIZATION EAM1 

1 PRINT 1000 EAM1 

CALL IRAN DPI 1»NMEAS» IA, IA, IA, IA, IA, I A, 4 ) EAM1 

CALL RANDPDI 1 , P SC AL E , D A , DA , D A , D A , D A , D A , A ) EAM1 

RETURN E AMI 

E AMI 

INTERROGATE FIGURE SENSOR AND SCALE MEASUREMENT E AMI 

SCALE POSITION EAM1 

2 DA = 0 . 0 E AMI 

DO 2000 1 = 1, NME A S EAM1 

X= X F SV ( I ) UPSCALE EAM1 

Y = YFSV( I UPSCALE EAM1 

E AMI 
E AMI 

INSERT DTOA AND ATOD SOFTWARE HERE EAM1 

( DA=DA+MEASURED FIGURE ERROR) EAM1 

EAM1 
E AMI 

2000 CONTINUE EAM1 

XFVI I )=DA*FSCALV( I ) /NMEAS E AMI 

RETURN EAM1 

C EAM1 

END E AMI 


1 5B0 
1590 
1600 
1610 
1620 
1630 
1 640 
1650 
1 660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 

17 60 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 

18 60 
1870 
1880 
1890 
1900 
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SUBROUTINE INI T ( NENTRY ) 

SUBROUTINE TO INITIALIZE THE EXPERIMENTAL ACTIVE MIRROR 

DIMENSION XVI 10000) , AMI 3000) ,AIM( 3000) , X FVI 20 ) , UF V ( 20 ) , A SC AL V ( 20 ) 

1 FSCALVI 20) ,XFSV< 20) , YFSVI 20) , XFRVI 20) ,DUMV< 20 ) , UF AV ( 20 ) , DUMVA ( 20 

2 »GA I NV( 10) ,GAINM< 400) , LACTV! 20) 

COMMON/BLKMFC/ N , NR , X V , MODE , L AC T V , MODOP , NSN SWT , NR A 
EQUIVALENCE ( XV( 1 ) ,AM( 1 ) ) , ( XVI 3001 ), AIM! 1 ) ) , ( XVI 6001 ) ,XFV( 1) ) , 

1 I UFVI 1 ) ,XV! 6101 ) ) ,< ASCALV! 1 > ,XV! 6121 )) ,( FSCALV! 1 ) ,XV< 6141 ) ) , 

2 I XFSVI 1) tXVI 6161) ) , I YFSVI 1) , XVI 6 1 8 1 ) ) * I XFRVI 1 ) , XV I 6201 ) ) , 

3 I UFA VI 1) tXVI 6221 ) ) ,< OUMVI 1 ) , XVI 6241 ) ) , I DUMVA I 1 ) , XV I 6261 ) > , 

4 I GAINVI 1 ) ,XV( 6281 ) ) ,( GAINMI 1 X XVI 6291 ) ) 

1000 FORMAT! 4HINIT) 

GO TO I 1,2) , NENTRY 

READ DATA FOR TYPEWRITER CONTROL 

1 PRINT 1000 

READ BASIC DATA FOR THE EXPERIMENTAL ACTIVE MIRROR 
N , NR 

CALL IRANDPI 2,N,NR*IA,IA,IA»IA,IA»4) 

AM 

CALL MXRNPI AM,N,N,4) 

FSCALV 

CALL MXRNPI FSCALV, 1,N, 4) 

XFSV 

CALL MXRNPI XFSV, 1,N, 4) 

YFSV 

CALL MXRNPI YFSV, 1,N, 4) 

LACTV 

CALL IMXRNPI LACTV, 1,N,4) 

A SC AL V 

CALL MXRNPI ASCALV, 1,NR,4) 

CALL FIGSENI 1, I ,XFV, XFSV, YFSV, FSCALV) 

CALL ACTCMDI 1 , I , UF V , UF AV , A SC ALV ) 

CALL MIRCALI 1 ) 

CALL ACTCALI 1 ) 

RETURN 

INITIALIZE MIRROR FIGURE CONTROL SYSTEM 

2 DO 2201 1=1, N 
X F V I I ) = 0.0 
XFRVI I )=0.0 
DUMVI I )=0.0 
DUMVA! I ) =0.0 

2201 UF V I I ) = 0.0 

DO 2202 1=1, NR 

2202 CALL AC TCMDI 2 , I , UF V , UF AV , A SC AL V ) 

DO 2203 1 = 1 , N 

2203 CALL F I GSENI 2 , I , XF V, XF SV , YF SV , F SC AL V ) 

RETURN 

END 


E AMI 19 1 0 
E AM 1 1 920 
E AMI 1930 
EAM1 1940 
, E AM 1 19 50 
) E AMI 1 960 
E AM 11970 
EAM11980 
EAM1 1990 
EAM12000 
EAM12010 
EAM12020 
E AM12030 
EAM12040 
EAM12050 
E AM 12 060 
EAM12070 
EAM12080 
EAM12090 
EAM12100 
EAM121 10 
EAM12120 
E AMI 2 1 30 
EAM12140 
EAM12150 
EAM12160 
EAM12170 
EAM12180 
E AMI 2 1 90 
EAM12200 
EAM12210 
EAM12220 
EAM12230 
E AMI 2240 
EAM12250 
EAM12260 
EAM12270 
EAM12280 
EAM12290 
E AMI 2 300 
EAM12310 
EAM12320 
EAM12330 
EAM12340 
EAM12350 
EAM12360 
EAM12370 
E AMI 2 380 
EAM12390 
EAM12400 
EAM12410 
EAM12420 
EAM12430 
EAM12440 
EAM12450 
EAM12460 
EAM12470 
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SUBROUTINE MFCS! NENTRY ) 

EAM12480 

c 


EAM12490 

c 

SUBROUTINE TO REALIZE MIRROR FIGURE CONTROL SYSTEMS 

EAM12500 

c 


EAM125 10 


DIMENSION XV! 10000) ♦ AM! 3000) , AIM! 3000) , XFV! 20 ) , UFV ( 20 ) , ASCALV ( 20 ) ,E AMI 2520 


1 FSCALV! 20) ,XFSV( 20) ,YFSVl 20) ,XFRV< 20) ,DUMV< 20) ,UFAV<20) ,DUMVAI 

1 20) EAM12530 


2 ,GAINV< 10) , GAINM! 400) , LAC TV! 20) 

EAM12540 


COMMON /RLK MFC/ N ,NR , X V , MODE » L AC TV , MODOP »NSN SWT t NR A 

EAM12550 


EQUIVALENCE ( X V! 1 ) , AM! 1 ) ) , ( XV! 3001) , AIM! 1 ) ) , ( X V ( 6001 ) , XFV ( 1 ) ) , 

EAM12560 


1 ( UFV! 1) ,XV( 6101 ) ) »( ASCALV! 1 ) , XV! 612 1 ) ) , ( FSCALV! 1 ) , XV ( 6141 ) ) , 

EAM12570 


2 ( XFSV! 1) ,XV< 6161) ) ,< YFSV! 1),XV< 6181) ) ,< XFRV! 1 ) ,XV< 6201 ) ) , 

EAM12580 


3 ( UFAV! 1 ) tXV! 6221 ) ) ♦ < DUMV! 1) ,XV( 6241 ) ) ,( OUMVA! 1 ) , XV ( 626 1 ) ) ♦ 

EAM12590 


4 ( GAINV! 1 ) ,XV( 6281 ) ) ,( GAINM! IX XV! 6291) ) 

EAM12600 


DIMENSION PARV! 20) ,GRAOV( 20) 

EAM12610 


COMMON/ BLKA/P INDEX, PI STOR , PARV, EPS» GR AOV, NI C, NHC , NP AR , N I M, NHM 

E AMI 2 620 

c 


EAM12630 


DIMENSION SPARV! 20 ) ,DV< 20) 

F AM 1 2 640 


COMMON /BLKA1 /NOPT ,EPSINC,EPSDEC,I SDEC»NSDEC, SNHC , SP AR V , DV 

EAM12650 

c 


EAM12660 

c 


EAM12670 

r 


EAM126B0 

1000 FORMAT! 4HMFCS) 

EAM12690 

1010 FORMAT! 1H1,/, 10X,5HGAINM) 

EAM12700 

c 


EAM12710 


GO TO ( 1, 2, 3, 4), NENTRY 

EAM12720 

C 


EAM12730 

C 

MIRROR FIGURE CONTROL SYSTEM INITIALIZATION 

E AMI 2740 

1 

PRINT 1000 

EAM12750 


CALL MXRNP! GAINV, 1,2,4) 

E AM 1 2 7 60 


CALL RANDPD! 1 , DT I TC S , D A , D A , 0 A , 0 A , DA , DA , 4 ) 

EAM12770 


I A = MODO P- 5 

EAM12780 

C 


E AM12790 

C 

CALCULATE FEEDBACK MATRIX 

EAM12ROO 


GO TO! 2110,2120,2130) , IA 

EAM12810 

c 


EAM12820 

c 

SLCS 

EAM12R30 

c 

GAINM=ARR**-1 

EAM12R40 

c 

GENERATE ARR 

EAM12850 

2110 NR A=NR 

E AM 1 2 8 60 


CALL REDUAM! 1 ) 

EAM12870 


CALL REDUAM! 2 ) 

EAM12880 

C 

GAINM=ARR**-1 

EAM12890 


CALL SINV! NR, AM, GAINM, DA ) 

E AMI 2900 

C 


EAM12910 


PRINT SLCS GAIN MATRIX 

EAM12920 


PRINT 1010 

EAM12930 


CALL MXRNP! GAINM, NR, NR, 3 ) 

EAM12940 


RETURN 

EAM12950 

c 


E AMI 2960 

C 

LOCS 

EAM12970 

C 

GA I NM = ART *AR 

E AM 1 2980 

2120 NR A=N 

EAM12990 


CALL MTRA! AM, AIM, N, NR, 0) 

FAM13000 


CALL MPRD! A I M , AM , GA I NM ,NR , N , 0, 0 ,NR ) 

EAM13010 

C 

GA I NM= ( ( ART*AR)**-1)*ART 

E AMI 3020 


CALL SINV! NR, GAINM, AM, DA) 

E AM13030 


CALL MPRD! AM, AIM, GAINM, NR, NR, 0,0, N) 

EAM13040 
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C 


EAM13050 


C PRINT LOCS GAIN MATRIX 

PRINT 1010 

CALL MXRNPi GA I NM , NR , N , 3 ) 

RETURN 

INITIALIZE ITERATIVE CONTROL SYSTEM 
2130 CONTINUE 
RETURN 

CALCULATE FIGURE CONTROLS 
ACTUATOR INPUT COMPUTATION 
MEASURE FIGURE ERROR 

2 DO 2505 1=1, N 

2505 CALL F I GSENI 2 , I , XFV , XF SV , YF SV , F SC AL V ) 

GENERATE XFRV 
I A=M0D0P-5 
J = 0 

DO 2506 1 = 1, N 

GO TO ( 2508,2509) ,1 A 

2508 IF< LACTVI I) 12507,2506,2507 
2507 J=J+1 

XFRV( J ) =XFV( I ) 

GO TO 2506 

2509 XFRV ( I ) = X F V ( I ) 

2506 CONTINUE 

DIJMV=GAINM*XFRV 

CALL MPRDI GAINM,XFRV, OUMV , NR, NR A, 0,0, 1 ) 

DUMV = DUMV*GAINV< 1 ) 

DO 2510 1=1, NR 

2510 DUMV( I )=DUMV< I )*GAINV( 1) 

INTEGRAL COMPENSATION 
DUMVA=DUMV*OT+DUMVA 

DO 2520 1=1, NR 

2520 DIJMVAI I )=DUMV( I )*DT+DUMVA( I ) 

ACTUATOR OUTPUT CONTROL SYSTEM 
UF V = GA I NV ( 2 ) *( DUMVA-UFAV !*DT-UFV 

3 DO 2530 1=1, NR 

2530 UF V( I ) =GA INV( 2)*( DUMVAt I )-UFAV( I ) )*DT + UFV( I ) 

TRANSFER COMMANDS TO ACTUATORS AND RETURN ACTUATOR OUTPUTS 
DO 2540 1=1, NR 

2540 CALL ACTCMDI 2 , I , UF V , UF AV , A SC ALV ) 

RETURN 

I TCS ACTUATOR COMMANDS 

ITCS COMMANDS 

4 CALL REALT( TREAL ) 

DA=TREAL 
DB=TREAL+DT ITCS 

C DUMVA=PARV 

DO 2550 1=1, NR 
2550 DUM V A ( I )=PARV( I ) 


EAM13060 
EAM13070 
E AM 1 3080 
EAM13090 
EAM13100 
EAM13110 
EAM13120 
E AMI 3 1 30 
E AMI 3 1 40 
FAM13150 
E AMI 3 1 60 
EAM13170 
EAM13180 
EAM13190 
EAM13200 
E AMI 32 10 
EAM13220 
EAM13230 
EAM13240 
EAM13250 
E AMI 32 60 
EAM13270 
E AMI 3280 
EAM13290 
EAM13300 
E AM 1 33 1 0 
EAM13320 
EAM13330 
EAM13340 
EAM13350 
E AMI 3 3 60 
EAM13370 
EAM13380 
E AM 1 3390 
E AMI 3400 
EAM13410 
E AMI 3420 
EAM13430 
E AM 1 3440 
EAM13450 
E AMI 3460 
EAM13470 
EAM13480 
E AM 1 3490 
EAM13500 
EAM135 10 
EAM13520 
EAM13530 
EAM13540 
EAM13550 
E AMI 35 60 
EAM13570 
EAM13580 
EAM13590 
EAM13600 
EAM13610 
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c 

C ACTUATOR OUTPUT CONTROL SYSTEM 

C UF V=GA I N( 2 ) *< DUMVA-UFAV > *OT-UFV 

2554 DO 2551 1=1, NR 

IJ F V ( I ) = GA INV( 2 ) ^( OUMV A( 1 )-UFAV( I ) )*DT + UFV( I ) 

2551 CALL ACTCMD( 2 , I , UFV , UF AV , A SC AL V ) 

DA=DA+DT 

2552 CALL REALT ( TREAL ) 

I F ( TREAL-OR >2555,2553, 2553 

2555 1 F ( TREAL-DA>2552, 2554, 2554 

2553 RETURN 
C 

END 


EAM13620 
EAM13630 
E AM 1 3 640 
EAM13650 
EAM13660 
EAM13670 
EAM13680 
EAM13690 
EAM13700 
E AMI 37 10 
EAM13720 
E AMI 3730 
EAM13740 
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SUBROUTINE MIR CAL I NENTRY ) 
SUBROUTINE TO CALIBRATE MIRROR 


EAM13750 
E AMI 37 60 
EAM13770 
EAM13780 


DIMENSION XVI 10000) , AMI 3000) , A I Ml 3000) ,XFV( 20) ,IJFV< 20) 

1 FSCALVI 20) , XFSVI 20) , YFSVI 20) ,XFRV< 20) , DUMVI 20) f UFAV<2 

2 tGAINVI 10)»GAINM( 400) ,L ACT VI 20) 

COMMON/BLKMFC/ N , NR , X V , MODE , L ACT V , MODOP , NSNSWT , NR A 
EQUIVALENCE (XVI 1 ) , AMI 1 > ) , I X VI 300 1 ) , A 1 Ml 1 ) ) ,( XVI 6001 ) , 

1 I UFVI 1 ) , XVI 6101 ) ) , I ASCALVI 1 ) ,XV< 6121 ) ) ,( FSCALVI 1 ) ,XV( 

2 I XFSVI 1) tXVI 6161) ) * I YFSVI 1) * X VI 6181) ) , I XFRVI 1) , XV I 620 

3 ( UFAVI 1 ) ,XV( 6221 ) ) ,< DUMMI 1 ) , XVI 6241 ) ) , I DUMVAI 1 ) ,XV( 62 

4 ( GAINVI 1) tXVI 62B1) ) , I GA1NMI 1 X XVI 6291 ) ) 


» ASCALVI 20) , EAM13790 
0) , DUMVAI 20 ) E AM 1 3800 
E AM 1 38 1 0 
EAM13820 
XFV(l)), EAM13830 

6141)), E AM 1 3840 

D), EAM13850 

61)), EAM13860 

FA M 13870 


E A M 1 3 8 8 0 


DIMENSION DUMMI 20) 

C 

1000 FORMAT! 6HMIRCAL ) 

1001 FORMAT! 1H1 / 1'3X , 2H AM ) 

1002 FORMAT! 8HMCAL END) 

C 

GO TO ( 1,2), NENTRY 
C 

C INPUT DATA 

1 PRINT 1000 

CALL IRANDPI 1 , NME AS , I A , I A , IA,IA,IA,IA,4) 

CALL RANDPDI 1 , DACT , DA , DA , DA , DA , DA , DA ,4 ) 

DB=2. 0*DACT 
RETURN 
C 

C MIRROR CALIBRATION 

2 CALL INITI2) 

DO 2200 1=1, NR 

DO 2201 J= 1 , NME AS 
UFVI I )=-DACT 

CALL ACTCMDI 2, I ,UFV,UFAV, ASCALV ) 

DO 2202 K= 1 , N 

CALL FIGSENI 2 , K , XF V , XF SV , YF SV , F SCAL V ) 

2202 DUMVI K )=XFV( K ) 

UFVI I ) = DACT 

CALL ACTCMDI 2 , I , UFV , UF AV , A SC AL V ) 

DO 2203 K= 1 , N 

CALL FIGSENI 2,K,XFV,XFSV,YFSV,FSCALV) 

DUMVI K ) = I X F V I K) -DUMVI K) ) /DB 

2203 AMI K+( I -1 ) *N ) = DUMVI K)+AM( K+( I-1)^N) 

2201 CONTINUE 

DO 2204 K= 1 , N 

2204 AMI K + l I — 1 )*N ) =AM( K + l 1-1 ) *N ) /NME AS 
2200 CONTINUE 

C 

C PRINT OUT MIRROR DEFORMATION-ACTUATOR COMMAND ARRAY 

PRINT 1001 

CALL MXRNPI AM, N, NR, 3) 

TYPE 1002 
RETURN 
C 

END 


EAM13890 
EAM13900 
E AM 1 39 1 0 
E AM 1 3920 
EAM13930 
EAM13940 
FAM13950 
E AMI 39 60 
EAM13970 
EAM13980 
E AMI 3990 
EAM14000 
E AM 1 40 1 0 
E AM 1 4020 
E AM14030 
E AMI 4040 
E AM 140 50 
E AMI 40 60 
E AMI 407 0 
EAM14080 
E AM14090 
EAM14100 
EAM141 10 
EAM14120 
EAM14130 
EAM14140 
EAM14150 
EAM14160 
EAM14170 
E AMI 4 180 
EAM14190 
E AMI 4200 
EAM142 10 
EAM14220 
EAM14230 
EAM14240 
EAM14250 
EAM14260 
E AM 1 42 70 
EAM14280 
EAM14290 
EAM14300 
EAM143 10 
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1000 

1010 

2020 

2010 

2030 

2040 


SUBROUTINE M00CHK! NENTRY , MODE , M0DED, M000P , M0D0PD, ITEST) 

SUBROUTINE TO CHECK VALUES OF MODE AND MODCIP 

FORMAT! 4HM0DE, 16, 5H N0T=,I5> 

FORMAT! 5HM0D0P, 15, 5H N0T=,I5) 

CHECK MODE 

IF! M0DE-M0DED)2020, 2010, 2020 
TYPE 1000, MODE, MOOED 
I TEST=2 

CHECK MODOP 

IF! MODOP-MODOPD) 2030,2040,2030 
TYPE 1010, MODOP, MODOPD 
I TE ST = 2 
RETURN 

END 


EAM14320 
E AMI 4330 
EAM14340 
EAM14350 
E AM 143 60 
EAM14370 
EAM143R0 
EAM14390 
EAM14400 
E AM 1 44 1 0 
EAM14420 
E AM14430 
E AM 1 4440 
EAM14450 
EAM14460 
EAM14470 
EAM14480 
E AM 1 4490 
EAM14500 


SUBROUTINE RE ALT ! TRE AL ) EAM14510 

EAM14520 

SUBROUTINE TO INTERROGATE REAL TIME CLOCK EAM14530 

EAM14540 

TREAL=RE AL TIME EAM14550 

EAM14560 

EAM14570 

INSERT REAL TIME CLOCK INTERROGATION SOFTWARE HERE EAM14580 

EAM14590 

EAM14600 

RETURN EAM14610 

EAM14620 

END EAM14630 
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SUBROUTINE REDUAMI NENTRY) 

SUBROUTINE TO GENERATE AR AND ARR FROM A 


E AM 1 4640 
EAM14650 
E AM 1 4660 
EAMI4670 


DIMENSION XV< 10000) ,AM(3000),AIM( 3000) ,XFV< 20), 

1 FSCALVI 20) ,XFSV< 20) ,YFSV( 20) »XFRV( 20) ,DUMV( 20) 

2 ,GAINV( 10) ,GAINM( 400),LACTV( 20) 

COMMON /BLKMFC/ N , NR , X V , MODE , L ACT V , M0D0P , NSN SWT , 
EQUIVALENCE I XVI 1 ) , AMI 1 ) ) , I XVI 3001 ) , AIM! 1 ) ) ,( XV 

1 (UFVI 1) ,XV( folOl ) )*(ASCALV( 1)»XV( 6121))»(FSCALV 

2 I XFSVI 1 ) ,XV( 6161 ) ) ,( YFSVI 1 ) ,XV( 6181 ) ) ,( XFRVI 1 ) 

3 I UF AVI 1 ) ,XV( 6 221 ) ) ,( DUMVI 1 ) ,XV( 6241 ) ) ,( DUMVAI 1 

4 ( GAINVI 1 ) ,XV( 6281 ) ) ,1 GAINMI 1 )( XVI 6291 ) ) 

1000 FORMAT! 1H 1 , / , 1 3 X , 2 HAR ) 

1010 FORMAT! 1H1 ,/, 12X,3HARR ) 

GO TOI 1,2) , NENTRY 

GENERATE AR BY REMOVING COLUMNS FROM A 
1 K = 0 

DO 2000 J=1,N 

IF( L ACT VI J) >2010,2000,2010 
2010 K = K+ 1 

DO 2020 1 = 1, N 

CALL LOCI I , J,IA,N,N,0) 

CALL LOCI I , K , IB,N,N,0) 


UFVI 20) , ASCALVI 20) ,E AM 14 680 
,UFAV( 20) , DUMVAI 20) E AMI 4690 

EAM14700 

NR A 

I 6001 ) , XFV I in. 


I 1) , XVI 6141 ) ) 
, XVI 6201 ) ) , 

) , X V I 6261 ) ) , 


EAM147 10 
E AM14720 
EAM14730 
EAM14740 
EAM14750 
E A M 1 4 7 60 
EAM14770 
EAM14780 
E AM 1 47 90 
EAM14800 
EAM14810 
E AM 14820 
EAM14830 
E AM 1 4840 
EAM14B50 
E AMI 48 60 
EAM14870 
E AMI 4880 
EAM14890 
EAM14900 


2020 AIM! I B ) = AM I IA) 


E AM 1 49 1 0 


2000 

CONTINUE 

EAM14920 

EAM14930 


COPY RESULT INTO AM 

E AM] 4940 


CALL MCPYI AIM, AM, N, NR, 0) 

EAM14950 

EAM14960 


PRINT AR 

E AM 1 4970 


PRINT 1000 

EAM14980 


CALL MXRNPI AM, N, NR, 3) 

E AM 1 4990 


RETURN 

E AMI 5000 
EAM15010 


GENERATE ARR FROM AR BY REMOVING ROWS FROM AR 

EAM15020 

2 

K = 0 

EAM15030 


DO 2100 1=1, N 

EAM15040 


IF! LACTVI I ) >2110,2100,2110 

E AMI 5050 

2110 

K = K+ 1 

E AMI 50 60 


DO 2120 J= 1 ,NR 

EAM15070 


CALL LOCI I,J, I A , N , NR , 0 ) 

EAM15080 


CALL LOCI K, J, IB, NR, NR, 0) 

EAM15090 

2120 

AIM! I B ) = AMI IA) 

EAM15100 

2100 

CONTINUE 

EAM151 10 
EAM15120 


COPY RESULT INTO AM 

EAM15130 


CALL MCPYI AIM, AM, NR, NR, 0) 

EAM15140 

EAM15150 


PRINT ARR 

E AMI 5 1 60 


PRINT 1010 

EAM15170 


CALL MXRNPI AM, NR, NR, 3) 

EAM15180 


RETURN 

EAM15190 

EAM15200 


END 

EAM15210 
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1000 

1001 

1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 
1011 
1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 


1 


2 


C 


SUBROUTINE TYPCON ! N6NTRY ) 

SUBROUTINE FOR TYPEWRITER CONTROL OF THE EXPERIMENTAL 

DIMENSION XVI 10000) , AMI 3000) , AIM! 3000) , XF V! 20 ) , UFV < 20 ) 

1 FSCALV! 20) ,XFSV( 20) ,YFSV( 20) ,XFRV< 20) ,DIJMV( 20) ,UFAV( 

2 tGAINVI 10)*GAINM( 400),LACTV< 20) 

COMMON /BLKMFC/ N , NR , XV , MODE , L ACTV , MODOP , NSN SWT , NR A 
EQUIVALENCE I XVI 1 ) ,AM( 1 ) ) , ( XVI 3001 ) , AIM! 1 ) ) ,< XVI 6001 ) ,XFV( 1) ) , 

1 ( UFVI 1 ) f XVI 6101 ) ) * ( ASCALVI 1 ) , XVI 612 1 ) ) , ( FSCALVI 1 ) ,XV I 6141 ) ) , 

2 ( XFSVI 1 ) ,XV{ 6161 ) ) ,( YFSVI 1 ) ,XV< 6181 ) ) , I XFRVI 1 ), XVI 6201 ) ) , 

3 ( UF AVI 1 ) ,XV( 6221 ) ) ,( DUMVI 1 ) ,XV! 6241) ) ,( DUMVAI 1 ) , XV I 6261 ) ) , 

4 I GAINVI 1 ) ,XV( 6281 ) ) , I GAINMI 1 ) I XVI 6291 ) ) 

DIMENSION PARVI 20) »GRADV( 20) 

COMMON/BLKA/PINDEX,PISTOR f PARV.EPS, GRADV,N I C,NHC,NPAR,NIM,NHM 
DIMENSION NAMVI 20) 

DATA NAMV / 4HXV ,4HAM ,4HAIM ,4HXFV ,4HUFV , 4HUF A V t 4H A SCV » 

1 4HFSCV,4HXFSV,4HYFSV,4HXFRV,4HDUMV,4HLACT,4HGANM,4HGANV,4HLACVt 

2 4HDEND/ 


FORMAT! A4 ) 

FORMAT! /6HNAMERR ) 
FORMAT! 1013 ) 

FORMAT! 3HNAM ) 
FORMAT! 4HINDX ) 
FORMAT! 10F12.6) 
FORMAT! 4HM0DE ) 
FORMAT! 4HINIT ) 
FORMAT! 4HSTRT ) 
FORMAT! 6HTYPC0N) 
FORMAT! 6HACTTST) 
FORMAT! 6HMIRTST ) 
FORMAT! 6HM0DFIN) 
FORMAT! 4HSLCS ) 
FORMAT! 4HL0CS) 
FORMAT! 4HITCS ) 
FORMAT! 1HJ) 

FORMAT! 2HJM) 

FORMAT! 4HCHNG ) 
FORMAT! 9HNEW VALUE ) 
FORMAT! 7HT00 BIG) 


GO TO ( 1 » 2 ) »NENTRY 

INITIALIZATION 
PRINT 1009 
NNAMV= 14 
I CHNG=2 
RETURN 

OPERATION 

CONTINUE 


EAM15220 
EAM15230 
ACTIVE MI RR0EAM15240 
EAM15250 
, ASCALV! 20) ,EAM15260 
20) *DUMVA( 20)EAM15270 
EAM15280 
EAM15290 
EAM15300 
EAM15310 
EAM15320 
EAM15330 
EAM15340 
EAM15350 
E AMI 5 360 
EAM15370 
FAM15380 
EAM15390 
E AMI 5400 
EAM15410 
EAM15420 
EAM15430 
EAM15440 
EAM15450 
E AMI 5460 
EAM15470 
EAM15480 
E A M 1 5 4 9 0 
E AMI 5500 
EAM15510 
EAM15520 
EAM15530 
EAM15540 
EAM15550 
EAM15560 
EAM15570 
EAM15580 
EAM15590 
EAM15600 
EAM15610 
EAM15620 
EAM15630 
EAM15640 
EAM15650 
EAM15660 
EAM15670 
EAM15680 
EAM15690 
EAM15700 
EAM157 10 
EAM15720 
EAM15730 
EAM15740 
EAM15750 
EAM15760 
EAM15770 
EAM15780 


74 



n n on on no on no noon on 



SELECT EXPERIMENTAL MODE 

EAM15790 

2280 

TYPE 1006 

EAM15800 


ACCEPT 1002, MODE 

EAM15810 


TYPE 1002, MODE 

EA M 15820 


GO TO (2211,2212,2213,2214,2215,2216,2217,2218,2219,2220,2221), 

EAM15830 

1 MODE 

EAM15840 

EAM15850 


INITIALIZE MFCS 

EAM15860 

2211 

TYPE 1007 

EAM15870 


GO TO 2207 

EAM15880 



EAM15890 


START MFCS 

EAM15900 

2212 

TYPE 1008 

EAM15910 


GO TO 2207 

EAM15920 

EAM15930 


DIAGNOSTIC MODE 

E AM 1 5940 

2213 

GO TO 2990 

EAM15950 
E AMI 5 9 60 


TEST ACTUATORS 

E AMI 5970 

2214 

TYPE 1010 

EAM15980 


CALL ACTCAL ( 2) 

EAM15990 


GO TO 2280 

EAM16000 
E AM 1 60 1 0 


TEST MIRROR 

EAM16020 

2215 

TYPE 1010 

EAM16030 


CALL MIRCALl 2) 

EAM16040 


GO TO 2280 

E AMI 6050 
EAM16060 


SIMPLIFIED LINEAR FIGURE CONTROL 

E AMI 6070 

2216 

TYPE 1013 

E AM 1 6080 


M0D0P=6 

EAM16090 


GO TO 2207 

EAM16100 
E AMI 6110 


LINEAR OPTIMAL FIGURE CONTROL 

EAM16120 

2217 

TYPE 1014 

EAM16130 


M0D0P=7 

E AMI 6140 


GO TO 2207 

EAM16150 
E AM 1 61 60 


ITERATIVE FIGURE CONTROL 

EAM16170 

2218 

TYPE 1015 

EAM16180 


M0D0P=8 

EAM16190 


GO TO 2207 

EAM16200 
E AMI 62 1 0 


EVALUATE AND TYPE FIGURE PERFORMANCE INDEX 

EAM16220 

2219 

CALL P I NDX ( 2 ) 

EAM16230 


TYPE 1016 

EAM16240 


TYPE 1005 , P I NDE X 

EAM1 6250 


DA=SORT( PINDEX/N) 

E AMI 62 60 


TYPE 1017 

EAM16270 


TYPE 1005, DA 

EAM16280 


GO TO 2280 

E AMI 6290 
EAM16300 


MODIFY DATA BUSS VALUE 

EAM16310 

2220 

ICHNG=2 

E AMI 6320 


TYPE 1018 

EAM16330 


GO TO 2990 

EAM16340 

EAM16350 
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C REQUEST MODE AGAIN IF MODE VALUE IS TOO LARGE 


EAM16360 


2221 TYPE 1020 
GO TO 2280 
C 

C IDENTIFY VARIABLE NAME 

2990 TYPE 1003 

ACCEPT 1000, KK 

TYPE 1000, KK 

DO 2302 1=1 , NNAMV 

I F( KK-NAMVl I ) >2302,2301,2302 

2301 LL= I 

GO TO 2303 

2302 CONTINUE 
TYPE 1001 
GO TO 2990 

2303 TYPE 1000, KK 

IDENTIFY VARIABLE INDEX 
TYPE 1004 

GO TO (2310,2320,2320,2310,2310,2310,2310,2310,2310,2310,2310, 
1 2310,2320,2310,2310,2280) ,LL 
2310 ACCEPT 1002,11 
TYPE 1002,11 
GO TO 2330 

2320 ACCEPT 1002, II, JJ 
TYPE 1002, II, JJ 
GO TO 2330 

ACCEPT NEW VALUE IF ICHNG=2 

2330 GO T0( 2331,2340) , ICHNG 

2331 ACCEPT 1005, V 
TYPE 1019 
TYPE 1005, V 
GO TO 2340 

TYPE VALUE OF INDEXED VARIABLE 

2340 GO TO (2401,2402,2403,2404,2405,2406,2407,2408,2409,2410,2411, 
1 2412,2413,2414,2415,2280) ,LL 

2401 CALL CHNG( ICHNG, V,XV( I I ) ) 

GO TO 2500 

2402 CALL ELMA( ICHNG, AM , I I , J J , V ,NR ) 

CALL ELMA( 2, AM, I I , J J , V ,NR ) 

TYPE 1005, V 

GO TO 2500 

2403 CALL ELMA< ICHNG, A IM, I I , J J, V,NR ) 

CALL ELMA( 2, AIM, I I ,JJ,V,NR) 

TYPE 1005, V 

GO TO 2500 

2404 CALL CHNG( ICHNG, V, XFV( I I ) ) 

TYPE 1005, XFV( I I ) 

GO TO 2500 

2405 CALL CHNG( ICHNG,V,XFRV( II)) 

TYPE 1005, XFRV( I I ) 

GO TO 2500 

2406 CALL CHNG( ICHNG, V , UFV( I I ) ) 

TYPE 1005, UFV( I I ) 

GO TO 2500 


EAM16370 
E AMI 6380 
E AMI 6390 
EAM16400 
E AMI 6410 
E AMI 6420 
E AMI 6430 
E AMI 6440 
EAM16450 
EAM16460 
EAM16470 
EAM16480 
E AMI 6490 
EAM16500 
EAM16510 
EAM16520 
E AMI 65 30 
E AMI 6540 
E AMI 6550 
E AMI 65 60 
EAM16570 
E AMI 6580 
E AMI 6590 
E AMI 6600 
E AM 1 661 0 
E AMI 6620 
EAM16630 
E AMI 6640 
EAM16650 
EAM1 6660 
E AM 1 667 0 
EAM16680 
E AMI 6690 
E AMI 6700 
EAM16710 
E AMI 6720 
EAM16730 
EAM16740 
E AMI 67 50 
E AMI 67 60 
EAM16770 
EAM16780 
E AMI 6790 
E AMI 6800 
EAM1 6810 
E AMI 6820 
EAM16830 
EAM16840 
EAM16850 
E AMI 68 60 
EAM16870 
E AMI 6880 
E AMI 6890 
EAM16900 
E AMI 69 10 
EAM16920 
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2407 

CALL CHNG( ICHNG*V,UFAV( II ) ) 

EAM16930 


TYPE 1005 ,UF AVt I I ) 

E AMI 6940 


GO TO 2500 

E AMI 6950 

2408 

CALL CHNG< ICHNG, V, ASCALVI I 1 ) ) 

EAM16960 


TYPE l/)05, ASCALVt I I ) 

E AM 1 6970 


GO TO 2500 

EAM16980 

2409 

CALL CHNGt I CHNG , V ,F SCAL V< 1 I)) 

EAM16990 


TYPE 1005,FSCALV< I I ) 

EAM17000 


GO TO 2500 

EAM17010 

2410 

CALL CHNG< ICHNG, V, XFSV< I I ) ) 

E AM 1 7020 


TYPE 1005 * XF SV( I I ) 

EAM17030 


GO TO 2500 

E AM 1 7040 

2411 

CALL CHNGt ICHNG, V,YFSV< II ) ) 

E AMI 7050 


TYPE 1005, YFSVt I 1 ) 

E AMI 70 60 


GO TO 2500 

EAM17070 

2412 

CALL CHNGt ICHNG, V, DUMVt I I ) > 

E AM 1 7080 


TYPE 1005, DUMVt I I ) 

EAM17090 


GO TO 2500 

EAM17100 

2413 

CALL ELMAt 2,GAINM, I I *JJ,V*NR) 

EAM17110 


TYPE 1005, V 

EAM17120 

2414 

TYPE 1 005 » GA I N V( II) 

EAM17130 


GO TO 2500 

E AM 1 7 140 

2415 

TYPE 1 002 , L ACT V( I I ) 

EAM17150 


GO TO 2500 

E AMI 7 1 60 

2500 

I CHNG= 1 

EAM17170 


GO TO 2280 

EAM17171 

2207 

RETURN 

EAM17180 

EAM17190 


END 

EAM17200 
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APPENDIX B 

PARAMETER OPTIMIZATION PACKAGE 

B. 0 Introduction 

This appendix contains listings for the following subroutines. 

Page 


SUBROUTINE E PCHNGI NENTRY , STEP ♦ STPOEC , D JOP V ) 80 
SUBROUTINE GRAD( NENTRY) 81 
SUBROUTINE I TPRTI NENTRY ) 82 
SUBROUTINE M INFCNl NENTRY ) 84 
SUBROUTINE P I NDXI NENTRY ) 86 

DUMMY SUBROUTINE PACKAGE 87 


SUBROUTINE ANGR ADI NENTRY > 
SUBROUTINE AVGR ADI NENTRY ) 
SUBROUTINE CNGRADI NENTRY > 
SUBROUTINE DA V I ONI NENTRY ) 
SUBROUTINE DESXI NENTRY) 
SUBROUTINE GRADMXI NENTRY ) 
SUBROUTINE M INF A( NENTRY ) 
SUBROUTINE NEWRAFI NENTRY ) 
SUBROUTINE PENFCNI NENTRY ) 
SUBROUTINE POWELI NENTRY ) 
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SUBROUTINE EPCHNG ! NENTR Y, STEP , S TPDEC , D JDP V ) 

SUBROUTINE TO LIMIT THE ABSOLUTE VALUE OF THE CHANGE IN THE 
COMPONENTS OF THE PARAMETER VECTOR 

DIMENSION D JDP VI 1 ) 

DIMENSION PARV! 20) ,GRADV( 20) 

COMMON /BLK A/P INDEX, PI STOR , PAR V, E P S, GR ADV, N I C , NHC, NPAR , N I M, NHM 
DIMENSION DPARLVI 20) , PARMINI 20) 

1000 FORMAT! 6HEPCHNG) 

GO TO! 1,2), NENTRY 

INPUT DATA 

1 PRINT 1000 

READ IN MAXIMUM ABSOLUTE CHANGE IN THE EACH ELEMENT OF THE 

PARAMETER VECTOR 

CALL MXRNP! DPARLV, 1, NPAR, 4) 

READ IN MINIMUM PARAMETER VALUE WITH SIGN 
CALL MXRNP! PARMIN, 1, NPAR, 4) 

RETURN 

SELECT EPS BASED ON ABSOLUTE VALUE OF CHANGE IN PARAMETER 

2 DO 2000 1=1 , NP AR 
2005 DA=DJDPV! I )*STEP 

IF! ABS! DA) -DPARLV! I ) 12020,2020,2010 
2010 NHC=NHC+1 

STEP=STEP*STPDEC 
GO TO 2005 
C 

2020 DB=P AR V! I ) -DA 

IF! P ARM IN! I ) 12021,2000,2 021 

2021 IF! ( DB-PARMIN! 1 ) )*PARMIN( I ) /ABS! PARMINI I ) ) ) 2030, 2000, 2000 
2030 DJDPV! I )=0.0 

PARV! I ) = PARM IN! I ) 

2000 CONTINUE 

CALL MFCS! 4 ) 

RETURN 

C 

END 


POPIOOOO 
POPIOOIO 
POP10020 
POP 10030 
P0P10040 
POP 10050 
P0P10060 
P0P10070 
POP 10080 
POP 1 0090 
P0P10100 
P0P10110 
P0P10120 
P0P10130 
P0P10140 
POP 101 50 
P0P10160 
P0P10170 
P0P101B0 
P0P10190 
P0P10200 
POP 102 10 
POP 10220 
P0P10230 
POP 10240 
P0P10250 
P0P10260 
POP 102 70 
P0P10280 
PQP10290 
P0P10300 
POP10310 
P0P10320 
POP10330 
POP 10340 
POP10350 
POP 10360 
P0P10370 
POP 10380 
P0P10390 
P0P10400 
P0P10410 
POP 10420 
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SUBROUTINE GRAD(NENTRY) 

POP 1 OA 30 

c 


POP 1 0 AAO 

c 

SUBROUTINE TO CALCULATE GRADIENT VECTOR 

POP 1 0A50 

c 

SPECIAL VERSION FOR OPTIMIZING THE PARAMETERS IN A DYNAMICAL 

SYSTEP0P10A60 

c 


POP10A70 


DIMENSION PARV( 20) ,GRADV< 20) 

P0P10A80 


COMMON /BL K A /PI NOE X, PI STOR , P AR V , EP S , GR ADV , N I C , NHC , NP AR ,NIM,NHM 

POP 10A90 

c 


P0P10500 

1001 

FORMAT ( AHGRAD) 

POP 105 10 

c 


P0P10520 


GO TO (1,2) »NENTRY 

P0P10530 

c 


P0P105A0 

1 

PRINT 1001 

POP 10 5 50 


CALL RANDPDI 1 , DPAR , DA , DA , OA , DA, DA , DA , A ) 

P0P10560 


CALL IRANDPt 1 , NGR AD , J J , J J , J J , J J , J J , J J , A ) 

POP1O570 


CALL ANGRADI 1 ) 

POP 1 0580 


RETURN 

P0P10590 

c 


POP 1 0 600 

c 

IF SENSE SWITCH 5 IS SET CALCULATE THE GRADIENT NUMERICALLY 

P0P10610 

2 

I F ( SENSE SWITCH 5)3,7 

P0P10620 

7 

GO T0( 3, A, 5 ) , NGR AO 

POP 10630 

C 


POP 1 0 fcAO 

C 

GRADIENT EVALUATION BY FINITE DIFFERENCES 

POP 10650 

3 

CALL P I NDX ( 2 ) 

POP 1 0 660 


DA = P I NDE X 

POP 10670 


00 1000 I = 1 , NPAR 

POP 1 0 680 


DB= P AR V( I ) 

POP 10890 


PARV( I )=DB+DPAR 

P0P10700 


CALL MFC S( A ) 

POP 1 07 1 0 


CALL P I NDX ( 2 ) 

P0P10720 


GR ADV ( I ) =( PINDEX-OA ) /DPAR 

POP 1 07 30 


PARV( I ) =DB 

P0P107A0 

1000 

CALL MFCS! A) 

P0P10750 


RETURN 

P0P10760 

C 


POP 10770 

C 

GRADIENT COMPUTATION BY TWO POINT INTERPOLATION 

P0P10780 

A 

DC=2.0*DPAR 

POP 1 0790 


DO 1100 1=1, NPAR 

POP 10800 


DB=PARV( I ) 

POP 108 10 


PARV( I ) =DB+DPAR 

P0P10820 


CALL MFCS! A) 

P0P10830 


CALL PINDXI2) 

POP 108A0 


DA=P I NDEX 

P0P10850 


PARV( I ) =DB-DP AR 

POP 1 08 60 


CALL MFC S( A ) 

P0P10870 


CALL PIN0X12) 

P0P10880 


GR ADV ( I )=( DA -P INDEX ) /DC 

P0P10890 


PARVl I )=DB 

P0P10900 

1100 

CALL MFCS! A) 

POP 10910 


RETURN 

POP 10920 

C 


P0P10930 

C 

CALCULATE GRADIENT ANALYTICALLY 

POP 1 09A0 

5 

CALL ANGRAD( 2) 

P0P10950 


RETURN 

POP 1 09 60 

C 


P0P10970 


END 

POP 10980 
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1000 

1010 

1040 

1050 

1060 


1 

2008 

2009 

2020 


2010 

2001 


2000 


110 


SUBROUTINE I TPRT < NENTRY ) 

SUBROUTINE TO PRINT OPTIMIZATION STATUS 
DIMENSION PARV! 20) ,GRADV( 20) 

COMMON/ BLKA/P INDEX, PI STOR , PARV , E P S , GRADV, N I C , NHC , NP AR , N I M , NHM 
DIMENSION SPI V( 20) ,SLGV( 20) 

FORMAT ( /, 12X,3HNIC, 1 2X , 3HNHC, 9X , 6HP I STOR, 9X , 6HP I NDEX , 9X , 6HDPI /PI 
I 12X,3HEPS,4X, 11H/GRADV/PARV) 

FORMAT ( 2115, 4F 15. 6) 

FORMAT! 5X»10H/SPIV/SLGV»6X,4HKKS=,I5! 

FORMAT! 3HSS3 ) 

FORMAT! 5HITPRT) 

GO TO < 1 ,2,3, 4,5) , NENTRY 

ITERATION OUTPUT 
I A = N I C + NHC 

IF((NIC-NIM)*( NHC-NHM) >2008,2020,2008 
IF! SENSE SWITCH 2)2020,2009 
IF! IA-IKS12010, 2020, 2020 
I KS= I A+NI TPRT 

PRINT N I C,NHC, PI STOR, PI NDEX, DP I /PI , EPS , GRADV , PARV 
PRINT 1000 

COMPUTE FRACTIONAL CHANGE IN PREVIOUS VALUE OF PERFORMANCE 
DA=( P INDEX-PI STOR ) /PI STOR 
PRINT 1Q10,NIC»NHC»PIST0R»PINDEX,DA,EPS 
CALL MXRNP(GRADV,1,NPAR,3) 

CALL MXRNP! P AR V , 1 ,NP AR , 3 ) 

CALL TO PINDX FOR AUXILLIARY OUTPUT PRINT 
CALL PINDX! 3 ) 

PRINT OUT MODIFIED PERFORMANCE INDEX WEIGHTING FACTORS 
CALL DESX! 3 ) 


INDEX 


KK = N I C+l 

IF! KK-KKS)2000, 2000, 2001 
I A=KKS-20 

STORED INFORMATION OUTPUTED EVERY 20 ITERATIONS 

PRINT 1040, IA 

CALL MXRNP! SPIV, 1 ,KKA, 3 ) 

CALL MXRNP! SLGV, 1, KKA,3 ) 

KKS=KKS+20 
KKA=KK— KKS+20 
SPIV! KKA ) =P I NDEX 

GENERATE THE LENGTH OF THE GRADIENT VECTOR AND STORE IT 
DA = 0. 0 

DO 2110 1 = 1 , NP AR 

DA = DA+GRADV! I )*GRADV( I ) 

SLGV! KKA ) = SORT < DA) 

GENERATE THE MAXIMUM ELEMENTS OF THE PARAMETER 
SET NIC=NIM+1 TO TERMINATE RUN IF SENSE SWITCH 
SET NHC=NHM+1 TO TERMINATE RUN IF SENSE SWITCH 
IF! SENSE SWITCH 3)2120,2300 


AND GRADIENT 
3 IS SET 
3 IS SET 


P0P10990 
POP110OO 
P0P11010 
P0P11020 
POP 1 1030 
P0P11040 
P0P11050 
POP11060 
P0P11070 
POP 11080 
POP 1 1090 
POPlllOO 
P0P11110 
P0P11 120 
POP 11130 
POP 1 1 140 
POP 1 1 1 50 
POP 1 1160 
POP 11170 
P0P11180 
POP 1 1190 
POP 1 1200 
POP 11210 
POP 1 1220 
P0P11230 
POP 1 1240 
POP 112 50 
POP 112 60 
P0P11270 
POP 11280 
POP 11290 
POP 1 1 300 
P0P11310 
POP 1 1 320 
P0P11330 
POP 1 1 340 
POP 1 1350 
POP 1 1360 
POP 11370 
POP 1 1 380 
P0P11390 
POP 1 1400 
POP 11410 
P0P11420 
P0P11430 
POP 1 1 440 
POP 1 1450 
POP 1 1460 
P0P11470 
POP 1 1 480 
POP 1 1490 
P0P11500 
POP 11510 
VECTORPOP1152® 
P0P11530 
POP 1 1540 
P0P11550 
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2120 TYPE 1050 


P0P1 1560 


PRINT 1050 
NIC=NIM+1 
NHC=NHM+ 1 
2300 RETURN 
C 

2 CONTINUE 
RETURN 

C 

3 CONTINUE 
RETURN 

C 

4 CONTINUE 
RETURN 

ITPRT INITIALIZATION 

5 CONTINUE 
PRINT 1060 

CALL IRANDPI 1,N ITPRT, I A, IA,IA,IA,IA,IA,4) 
I KS = 0 
KKS=20 
RETURN 
C 

END 


P0P1 1570 
POP 1 1 580 
P0P1 1590 
P0P11 600 
POP 1 1610 
P0P1 1 620 
POP 1 1 630 
P0P11 640 
P0P1 1650 
P0P11 660 
P0P1 1670 
P0P1 1 680 
POP 1 1 690 
POP 1 1700 
POP 1 17 10 
POP 11720 
POP 11730 
POP 1 1 740 
P0P11750 
POP 1 1 7 60 
P0P11770 
POP11730 
P0P1 1790 
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SUBROUTINE MINFCN! NENTRY) 


P0P11890 


1000 

1001 

1120 


11 


12 


13 


C 

C 


900 

C 

C 


SUBROUTINE TO MINIMIZE A PERFORMANCE INDEX 

SPECIAL VERSION FOR OPTIMIZING THE PARAMETERS IN A DYNAMICAL 
DIMENSION PARV! 20) ,GRADV< 20) 

COMMON/BLKA/P INDEX, PI ST0R,PARV, EPS,GRADV,NIC,NHC,NPAR,NIM,NHM 
DIMENSION SPARV! 20) ,DV( 20) 

C0MM0N/BLKA1/N0PT,EPSINC,EPSDEC,I SDEC,NSDEC,SNHC,SPARV,DV 

FORMAT! 7E10.0) 

FORMAT! 7F15.6) 

FORMAT! 6HMINFCN) 

GO TO! 11,12,13) , NENTRY 

INPUT DATA 
PRINT 1120 
CALL RANDPI2) 

CALL RANDPD! 3 , E PS, EPS INC , EP SDEC, DA , DA , DA, DA, 4 ) 

CALL IRANDP! 6 , N I M ,NHM,NOP T , NPAR , I SDEC,NSDEC, JJ,4) 

EPSM=EPS 


P0P1181G 
POP 11820 
SY STEPOP 1 18 30 
POP 1 1 840 
P0P11850 
POP 118 60 
POP 11870 
POP 1 1 880 
P0P11890 
POP 1 1 900 
POP 1 1910 
POP 1 1920 
POP 1 1930 
POP 1 1 990 
POP 1 1950 
POP 1 1 960 
POP 11970 
POP 1 1 98© 
POP 1 1990 
P0P12000 
POP1201G 
POP 12020 


OPTIMIZATION ALGORITHMS 
PERFORMANCE INDEX INITIALIZATION 
CALL IRANDP! 1 , NP AR , I A , I A , I A , I A , I A , 1 A, 4 ) 
CALL MXRNP! PARV, 1, NPAR, 4) 

CALL GRAD! 1 ) 

CALL EPCHNG! 1 , EPS, EPSDEC,GRAOV ) 

CALL ITPRTI5) 

CALL GRADMX! 1 ) 

CALL MINFA! 1) 

CALL AVGRAD! 1 ) 

CALL NEWRAF! 1) 

CALL CNGRAD! 1 ) 

CALL DAVIDN! 1) 

CALL POWEL! 1 ) 

RETURN 

INITIALIZE PERFORMANCE INDEX 
CALL PINDX! 1 ) 

RETURN 

PERFORM MINIMIZATION 

PIST0R=10.0**6 

N I C = — 1 

NHC = 0 

SNHC=-1 

EPS=EPSM 

INITIALIZE UFV 

CALL MFCS! 9 ) 

INITIALIZE GRADV , GRAOM 
DO 900 1=1, NPAR 
GRADV! I ) =0 .0 

GO TO! 1,2, 3, 4, 5, 6) ,NOPT 


POP12030 
POP12040 
POP1205O 
P0P12Q6Q 
P0P12070 
POP12080 
P0P12090 
POP 1 2 100 
POP 12110 
P0P12120 
P0P12130 
P0P12 140 
P0P12150 
POP 12 1 60 
P0P12170 
POP 12 180 
P0P12190 
P0P12200 
P0P12210 
POP 12220 
P0P12230 
P0P12240 
P0P12250 
P0P12260 
P0P12270 
POP 1 2280 
P0P1229O 
P0P12300 
P0P12310 
P0P12320 
P0P12339 
P0P12340 
P0P12350 
P0P12360 
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c 

METHOD OF STEEPEST DESCENT 

POP 12 370 

1 

GO TO 1105 

POP 1 2 380 

c 

BEGINNING OF ITERATIVE LOOP 

POP 1 2 390 

c 

STORE PINDEX AND PARV 

P0P12400 

1011 

P I STOR = P INDEX 

P0P12410 


DO 1002 1 = 1* NP AR 

P0P12420 

1002 

SPARVI I )=PARV( I ) 

P0P12430 

C 

AUXILLIARY STORAGE 

P0P12440 


CALL P I N DX ( 4 ) 

PGP 1 2450 


CALL GR AD( 2) 

P0P12460 

c 

MODIFIED STEEPEST DESCENT 

P0P12A70 

c 

SEARCH FOR A MINIMUM ALONG THE DIRECTION OF THE GRADIENT VECTOR 

POP 1 2 ARO 

c 

IF NSDEC = 2 

P0P12490 


GO T0( 1115,1121 ) *NSDEC 

P0P12500 

1121 

CALL M INF A( 2) 

P0P125 10 


GO TO 1105 

P0P12520 

1115 

CALL EPCHNGI 2 * E P S ,E P SDEC * GR ADV ) 

P0P12530 


DO 1100 I=1*NPAR 

P0P12590 

1100 

PARV( 1 ) = SP AR V( I )-EPS*GRADV( I ) 

P0P12550 


CALL MFCS(A) 

POP 125 60 

1105 

CALL PINDXI2) 

P0P12570 


1 F ( P I NDEX-PISTOR) 1101, 1102, 1102 

P0P12580 

1101 

N I C=N I C+ 1 

P0P12590 


I F C NHC-SNHC) 1200, 1200, 1210 

POP 1 2 600 

1200 

EPS=EPS*EPSINC 

POP 12 610 

1210 

SNHC=NHC 

P0P12620 


CALL I TPRT ( 1) 

POP 1 2 630 


I F ( NIC-NIM) 1011, 1110,1110 

POP 1 2 640 

1102 

NHC=NHC+1 

POP 1 2 650 


EPS=EPS*EPSDEC 

POP 1 2 660 


CALL I TPRT ( 1) 

POP 1 2 67 0 


I F( NHC-NHM) 1115, 1110, 1110 

POP 1 2 680 

1110 

RETURN 

POP 1 2 690 

C 


P0P12700 

C 

METHOD OF AVERAGE GRADIENT 

P0P127 10 

2 

CONTINUE 

P0P12720 


CALL AVGRADI 2 ) 

P0P12730 


RETURN 

POP 12740 

C 


POP12750 

C 

METHOD OF NEWTON RAPHSON 

POP 1 27 60 

3 

CONTINUE 

POP12770 


CALL NEWRAFI 2 ) 

P0P12780 


RETURN 

POP 1 2 790 

C 


POP 1 2 800 

C 

METHOD OF THE CONJUGATE GRADIENT 

P0P128 10 

4 

CONTINUE 

P0P12820 


CALL CNGRADt 2) 

P0P12830 


RETURN 

POP 1 2840 

C 


POP 1 2850 

C 

METHOD OF DAVIDON 

POP 12860 

5 

CONTINUE 

P0P12870 


CALL DA V I DN( 2 ) 

POP128RO 


RETURN 

P0P12890 

C 


POP12900 

C 

METHOD OF POWELL 

P0P12910 

6 

CONTINUE 

P0P12920 


CALL POWEL ( 2 ) 

POP12930 


RETURN 

POP 12940 

C 


POP 12950 


END 

POP 12960 
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on non noon 


SUBROUTINE PINDXI NENTRY ) 


POP 12970 


P0PI2980 

SUBROUTINE TO EVALUATE FIGURE PERFORMANCE P0P12990 

INDEX J = XT*X POP 1 3000 

POP 1 30 1 0 

DIMENSION XV{ 10000) ,AM< 3000) , AIM( 3000) ,XFV(50),UFV(50),ASCALV<50> , POP 13020 
1 FSCALVI 50) , XFSV( 50 > ,YFSV( 50) ,XFRV< 50) ,DUMV( 50) ,UFAV< 50) , L ACTV ( 50 ) POP 1 3030 
COMMON/BLKMFC/ N , NR , NR A, XV , MODE , L ACTV , MODOP , NSNSWT POP13G40 

EQUIVALENCE ( XV( 1 ) ,AM( 1 ) ) ,{ XVI 3001 ) ,AIM( 1 ) ),( XVI 6001 ) ,XFV( 1 ) ) , POP 13050 

1 ( UFV( 1 ) ,XV( 6101) ) ,( ASCALVI 1 ) , X V( 61 5 1 ) ) » ( FSCALVI 1 ) ,XV( 6201 ) ) , POP 13060 

2 < XFSV(l) ,XV< 6251 ) ) , ( YFSV< 1 ) ,XV( 6301 ) ) ,( XFRVt 1 ) ,XV( 6351 ) ) , POP 13070 

3 ( UFAVI 1) ,XV< 6401) ),( DUMV( 1 ),XV( 6451) ) P0P13080 

POP 1 3090 

GO T0( 1 1 2 ) tNENTRY P0P13100 

POP 13110 


INITIALIZATION 
1 CONTINUE 


P0P13120 
POP 13130 


RETURN 


POP131AO 


P0P13150 

CALCULATE PINDEX=XFVT*XFV P0P13160 

2 PINDEX=0.0 P0P13170 

DO 2100 1=1, N P0P13180 

CALL FIGSENI 2 , I , XFV , XFSV , YF SV , F SCAL V ) P0P13190 

2100 PINDEX=PINDEX+XFV( I)*XFV< I) P0P13200 

RETURN POP132 10 

C POP 1 3220 

END P0P13230 
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DUMMY SUBROUTINE PACKAGE 


SUBROUTINE 

ANGRADI NENTRY ) 

P0P13240 

RETURN 


POP 132 50 

END 


POP 1 3 2 60 

SUBROUTINE 

AVGRADI NENTRY ) 

P0P13270 
POP 1 32 80 
POP I 3290 

RETURN 
END ' 


SUBROUTINE 

RETURN 

END 

CNGRADI NENTRY ) 

P0P13300 
POP 1 33 1 0 
POP 13320 

SUBROUTINE 

DAVIONI NENTRY ) 

P0P13330 

RETURN 

POP 13340 

END 


POPI 3350 

SUBROUTINE 

DESX( NENTRY) 

POP 13360 

RETURN 


P0P13370 

END 


POP 1 3380 

SUBROUTINE 

GRADMXI NENTRY ) 

POPI 3390 

RETURN 


POPI 3400 

END 


POP 13410 

SUBROUTINE 

MINFAt NENTRY) 

P0P13420 

RETURN 


P0P13430 

END 


POPI 3440 

SUBROUTINE 

NEWRAFI NENTRY ) 

POP 13450 

RETURN 


POP 1 3460 

END 


POP 1 3470 

SUBROUTINE 

PENFCNI NENTRY ) 

P0P134B0 
POP 13490 

RETURN 

END 


POP 1 3500 

SUBROUTINE 

POWEK NENTRY ) 

POP 1 35 1 0 

RETURN 


POP 13520 

END 


POP 1 35 30 
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preceding page 


blank not filmed 


APPENDIX C 

MATHEMATICAL OPERATIONS PACKAGE 


C. 0 Introduction 


This appendix contains listings for the following programs. 


Page 


FUNCTION ELM(A,L,M,N> 90 
SUBROUTINE ELMA ( NENTR V , A , I , J , V , N ) 90 
SUBROUTINE GMADDI A , B , R , N , M ) 91 
SUBROUTINE GMPROI A , B , R , N, M , L ) 91 
SUBROUTINE GMSUBI A , B , R , N , M ) 92 
SUBROUTINE GMTR A( A , R , N» M ) 92 
SBUROUT I NE GTOSYMI X , X S ,NX ) 93 
SUBROUTINE LOCI I , J, IR,N,M,MS) 93 
SUBROUTINE MCP Y( A , R , N , M , MS ) 94 
SUBROUTINE MMADDI N, ALPHA, AtBETA, B,C ) 94 
SUBROUTINE MPRUI A , B ,R , N , M ,M S A , MSB , L ) 95 
SUBROUTINE MTRA I A , R , N , M , MS ) 96 
SUBROUTINE SYMTOGI X S , X ,NX ) 96 
SUBROUTINE S I NVI N , A I , B, 0 > 97 
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FUNCTION ELM ( A,L,M,N) 

FUNCTION RETURNS THE VALUE OF THE L * M TH ELEMENT OF THE MATRIX 
A WHICH HAS N COLUMNS AND AN ARBITRARY NUMBER OF ROWS 

DIMENSION A( 1 ) 

ELM=A( M*N-N+L ) 

RETURN 

END 


MOP 10000 
MOP 1 001 0 
MOP 1 0020 
M0P10030 
M0P100A0 
MOP 10050 
MOP 10060 
M0P10070 
M0P10080 
MOP 1 0090 
MOP 101 00 


1 


2 

C 

c 


SUBROUTINE ELMA( NENTRY.A. I . J.V.N) 

SUBROUTINE TO WRITE INTO ANO READ FROM MEMORY THE I,JTH ELEMENT 
OF MATRIX A WHICH IS STORED IN GENERAL FORM. 

DIMENSION A ( 1 ) 

GO TOI 1,2) ,NENTRY 

A ( I,J) = V 

A ( I+( J-l )*N)=V 

RETURN 

V= A ( I , J ) 

V = A ( I + ( J-l )*N) 

RETURN 

END 


M0P10110 
M0P10120 
M0P10130 
MOP 101 AO 
M0P10150 
MOP 10160 
MOP 10 1 7 0 
MOP 101 BO 
MOP 10190 
M0P10200 
M0P10210 
MOP 10220 
M0P10230 
MOP 10240 
MOP 10250 
MOP 102 60 
M0P10270 
MOP10280 
MOP 1 0290 
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SUBROUTINE GMADD ( A,B.R.N.M) 


SUBROUTINE PERFORMS MATRIX ADDITION, R=A+R» WHERE A , R AND R ARE 
N BY M MATRICES. 

DIMENSION A( 1 ) , B( 1 ) , R( 1 ) 

NM=N*M 

DO 110 1 = 1, NM 
110 R ( I ) = A ( I ) +B( I ) 

RETURN 

END 


MOP 1 0300 
MOP10310 
M0P10320 
MOP 10330 
MOP 10340 
M0P10350 
MOP 1 03 00 
MOP 1 0370 
MOP 103BO 
MOP 1 0 390 
MOP 1 0400 
MOP 1 0 A 1 0 
MOP10420 


SUBROUTINE GMPRD ( A,B,R,N,M,L) M0P10430 

M0P10440 

FORM THE PRODUCT R=A*B WHERE A IS A N*M MATRIX AND R IS A M*L MATRM0P1 0450 
A , B AND R ARE STORED IN GENERAL MATRIX FORM COLUMN BY COLUMN M0P10460 


MOP 10470 


DIMENSION A( 1) ,B( 1) ,R( 1) 

MOP 10480 


MOP 1 0490 

I R = 0 

MOP 1 05 00 

I K = — M 

MOP 1 05 1 0 

DO 10 K= 1 , L 

MOP 10520 

I K= I K+M 

M0P10530 

DO 10 J= 1 , N 

MOP 10540 

I R= I R+ 1 

MOP 1 0 550 

J I = J-N 

MOP 105 60 

I B= I K 

MOP 10570 

R ( IR) =0.0 

MOP 10580 

DO 10 I = 1 , M 

MOP 1 0590 

J I = J I +N 

MOP 10600 

I fi= I B+ 1 

M0P10610 

R( I R ) =R( IR ) +A( J I ) *B( IB) 

MOP 1 0620 

RETURN 

MOP 1 0 630 


MOP 1 0 640 

END 

MOP 1 0 650 
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SUBROUTINE GMSUBl A , R , R , N , M ) 

SUBROUTINE PERFORMS MATRIX SUBTRACTION, R = A-B , 
N BY M MATRICES. 

DIMENSION A( 1 ) ,R( 1 ) ,R( 1 ) 

NM = N*M 

00 110 1=1, NM 
110 R( I) = A( I ) — B ( I) 

RETURN 

ENO 


WHERE A , B AND 


MOP 10660 
MOP 1 0 670 
R AREM0P10680 
MOP 1 0 690 
MOP 10700 
MOP 1 07 1 0 
MOP 107 2 0 
MOP 1 0730 
MOP 10740 
M0P10750 
MOP 1 07 60 
MOP 10770 
MOP 107 R 0 


SUBROUTINE GMTRA ( A , R , N , M ) 

TRANSPOSE A GENERAL MATRIX 

A - NAME OF MATRIX TO BE TRANSPOSED 
R - NAME OF RESULTANT MATRIX 
N - NUMBER OF ROWS OF A AND COLUMNS OF R 
M - NUMBER OF COLUMNS OF A AND ROWS OF R 

DIMENSION A( 1) ,R( 1) 

C 

I R = 0 

DO 10 1=1, N 
I J= I -N 
DO 10 J= 1 , M 
I J=I J + N 
I R= I R+ 1 

10 . R( I R ) =AI I J ) 

RETURN 

C 

END 


MOP 1 0790 
MOP 1 OB 00 
M0P10810 
M0P10820 
MOP 10830 
MOP 10840 
MOP 10850 
M0P10860 
M0P10870 
M0P10880 
MOP 10890 
MOP 1 0900 
MOP 1 09 1 0 
M0P10920 
MOP 1 0930 
MOP 10940 
MOP 1 0950 
MOP 109 60 
M0P10970 
MOP 10980 
MOP 1 0990 
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10 


SUBROUTINE GTOSYMI X , X S , NX ) 


PROGRAM CON V 
NF * ( NF+1 ) /2 
THE NX BY NX 
THE NX BY NX 
IN COLUMNAR 


ERTS A SQUARE NX BY NX MATRIX INTO A VECTOR WHOSE 
AND WHOSE ELEMENTS CONSIST OF THE UPPER TRIANGLE 
MATRIX? STORED IN COLUMNAR FORM. 

MATRIX MUST BE STORED IN A VECTOR WHOSE LENGTH I 
FORM, BEFORE THIS ROUTINE IS CALLED. 


DIMENSION X( 1 ) ,XS( 1 ) 
L L = 0 

DO 10 J = 1 , N X 
DO 10 I=1,J 
LL=LL+1 
K=( J-l )*NX+I 
XS( L L ) = X( K ) 

RETURN 


END 


MOP 

1 

1000 

MOP 

1 

1010 

LENGMOP 

1 

1020 

OF MOP 

1 

1030 

MOP 

1 

1040 

S NX *MOP 

1 

1050 

MOP 

1 

1060 

MOP 

1 

1070 

MOP 

1 

lOBO 

MOP 

1 

1090 

MOP 

1 

1100 

MOP 

1 

1110 

MOP 

1 

1120 

MOP 

1 

1130 

MOP 

1 

1140 

MOP 

1 

1150 

MOP 

1 

1160 

MOP 

1 

1 170 

MOP 

1 
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SUBROUTINE LOC( I , J , I R , N , M , MS ) 

MOP 11190 
MOP 11200 


SUBROUTINE TO GENERATE VECTOR SUBSCRIPT FOR AN ELEMENT IN A 

MATRIXM0P11210 


OF SPECIFIED STORAGE MODE. 

M0P1 1220 


MS =0 SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N*M ELEMENTS 

MOP 11230 


IN STORAGE (GENERAL MATRIX) 

MOP 1 1240 


MS=1 SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N*(N+l)/2 IN 

MOP 1 1 250 


STORAGE (UPPER TRIANGLE OF SUMMETRIC MATRIX). IF 

M0P1 1260 


ELEMENT IS IN LOWER TRIANGULAR PORTION, SUBSCRIPT IS 

MOP 11270 


CORRESPONDING ELEMENT IN UPPER TRIANGLE. 

M0P1 1280 


MS=2 SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N ELEMENTS 

MOP 1 1 290 


IN STORAGE (DIAGONAL ELEMENTS OF DIAGONAL MATRIX). 

MOP11300 


IF ELEMENT IS NOT ON DIAGONAL (AND THEREFORE NOT IN 

MOP 11310 


STORAGE), IR IS SET TO ZEROi 

MOP11320 
M0P1 1330 


IX=I 

MOP 1 1340 


J X = J 

MOP 1 1 350 


I A = I — J 

M0P11360 


IF(MS-l) 10,20,30 

M0P11370 

10 

IRX=N*( JX-1 )+IX 

M0P1 1380 


GO TO 36 

M0P1 1390 

20 

I F ( IA)22,24,24 

M0P1 1400 

22 

I RX= I X+ ( JX# JX-JX ) / 2 

MOP 1 1410 


GO TO 36 

MOP 1 1420 

24 

I R X = J X + ( I X* I X — I X )/ 2 

MOP 1 1 430 


GO TO 36 

M0P1 1440 

30 

I RX = 0 

M0P1 1450 


I F ( IX-JX >36,32,36 

MOP 1 1460 

32 

I RX= I X 

MOP 1 1 470 

36 

I R = I R X 

MOP 11480 


RETURN 

MOP 1 1 490 
M0P1 1500 


END 

M0P11510 
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1 


SUBROUTINE MCPYI A ,R , N , M , M S ) 


MOPI 1520 


MCPY COPIES ENTIRE N BY M MATRIX A INTO N BY M MATRIX R 

MS - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A (AND R) 

0 - GENERAL 

1 - SYMMETRIC 

2 - DIAGONAL 

DIMENSION A{ 1 ) , R ( 1 ) 

COMPUTE VECTOR LENGTH, IT 
CALL LOC( N,M, IT ,N,M»MS) 

COPY MATRIX 
DO 1 1=1, IT 
R( I ) =A( I ) 

RETURN 

END 


MOP 115 30 
MOP 1 1 5 AO 
MOPI 1550 
MOP 115 60 
MOP 11570 
MOPI 1580 
MOP 11590 
MOP 1 1 600 
MOP 11610 
MOP 1 1 620 
MOP 1 1 630 
MOP 1 1 640 
MOPI 1 650 
MOP 1 1 660 
M0P11 670 
MOP 11680 
MOP 1 1 690 


SUBROUTINE MM AOD ( N , ALPH A , A , BE T A , B , C ) 

SUBROUTINE TO FORM COMBINATION C=ALPHA*A+BETA*B 

DIMENSION A( 1) ,B( 1) ,C< 1) 

C 

DO 1 1=1, N 

1 C( I )' = ALPHA*A< I )+BETA*R( I ) 

RETURN 

C 

END 


MOP 117 00 
MOPI 1710 
MOPI 1720 
MOPI 1730 
MOP 11740 
MOP 11750 
MOP 117 60 
MOPI 1770 
MOP 117 80 
MOP 11790 
MOPI 1800 
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10 

20 


30 


40 


50 

60 

70 

80 

90 


SUBROUTINE MPRD ( A,B,R,N,M,MSA,MSR,L) 

MPRD MULTIPLIES N BY M MATRIX A BY M BY L MATRIX 
PRODUCT INTO N BY L MATRIX R 

MSA - ONE DIGIT NUMBER FOR STORAGE MODE OF 

0 - GENERAL 

1 - SYMMETRIC 

2 - DIAGONAL 

MSB - SAME AS MSA EXCEPT FOR MATRIX R 

DIMENSION A< 1 ) ,B( 1 ) f R( 1 ) 

SPECIAL CASE FOR DIAGONAL BY DIAGONAL 

MS=MSA*10.+MSR 

I F( MS-22 ) 30,10,30 

DO 20 1=1, N 

R< I ) = A < I )*B( I ) 

RETURN 


MOP 11810 

B AND STORES THE M0P11820 
M0P11830 

MATRIX A MOP 1 1840 

M0P1 1850 
M0P1 I860 
MOP 11870 
MOP 1 1880 
M0P1 1890 
M0P1 1900 
MOP 11910 
MOP 119 20 
M0P1 1930 
M0P1 1940 
M0P1 1950 
M0P1 1960 
MOP 1 1970 


ALL OTHER CASES 
I R= 1 

DO 90 K= 1 , L 

DO 90 J= 1 , N 

R ( I R ) = 0 

DO 80 1 = 1, M 

I F ( MS >40,60,40 

CALL LOCI J,I,IA,N,M, MSA ) 

CALL LOCI I , K , I B , M , L , M SB ) 

IF! I A ) 50 , 80 , 50 

IF( 18)70,80,70 

IA = N*< I-D+J 

IB=M*< K- 1 ) + I 

R I I R ) =R( IRl+AI I A ) *8 1 IB) 

CONTINUE 

I R= I R+ 1 

RETURN 

END 


MOP11980 
M0P1 1990 
MOP 12 000 
MOP 12 010 
MOP 1 2020 
MOP 1 2 030 
MOP 12040 
M0P12050 
MOP 12060 
MOP 12070 
MOP 12080 
MOP 12 090 
MOP 12100 
MOP 12110 
MOP 12120 
M0P12130 
MOP 1 2 1 40 
MOP 1 2150 
M0P121 60 
MOP 1 2170 
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SUBROUTINE MTR &( A ,R , N , M , MS > 

MOP 1 2 1 BO 

c 


MOP 12190 

c 

MTR A TRANSPOSES N BY M MATRIX A TO FORM M BY N MATRIX R 

M0P12200 

c 

MS - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A 

( AND R) M0P12210 

c 

0 - GENERAL 

M0P12220 

c 

1 - SYMMETRIC 

M0P12230 

c 

2 - DIAGONAL 

M0P12240 

c 


M0P12250 


DIMENSION A ( 1 ) * R ( 1 ) 

M0P12260 

c 


M0P12270 

c 

IF MS IS 1 OR 2, COPY A 

M0P122R0 


I F ( MS) 10,20,10 

M0P12290 

10 

CALL MCPY( A ,R ,N,N ,MS ) 

MOP 1 2300 


RETURN 

M0P12310 

c 


MOP 12 320 

c 

TRANSPOSE GENERAL MATRIX 

M0P12330 

20 

I R = 0 

MOP 12 340 


DO 30 1=1, N 

MOP 1 2 350 


I J= I -N 

M0P12360 


DO 30 J= 1 , M 

M0P12370 


I J=I J+N 

M0P12380 


I R= I R+ 1 

MOP 1 2390 

30 

R ( 1R ) =A( I J) 

MOP 12400 


RETURN 

M0P12410 

C 


M0P12420 


END 

M0P12430 



SUBROUTINE S YMTOGI X S , X , NX ) 

MOP 1 2440 

C 


M0P12450 

C 

PROGRAM CONVERTS A SYM. MATRIX VECTOR (IN SUPPRESSED SYM. 

STOR AGE ) MOP 1 2460 

C 

WHOSE LENGTH IS NX*(NX+l)/2, INTO A GENERAL MATRIX VECTOR 

WHOSE MOP 12470 

C 

LENGTH IS NX*NX. 

M0P124B0 

C 


MOP 1 2490 


DIMENSION X( 1) ,XS< 1) 

MOP 12500 

C 


M0P12510 


LL = 0 

M0P12520 


DO 10 J= 1 , NX 

M0P12530 


DO 10 1=1, J 

M0P12540 


LL=LL+1 

M0P12550 


K=( J-l )*NX+I 

M0P12560 


M= ( I — 1 ) *NX + J 

M0P12570 


X( M)=XS( LL ) 

M0P12580 

10 

X( K ) = X S ( LL ) 

M0P12590 


RETURN 

MOP 12 600 

C 


MOP 12 610 


END 
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SUBROUTINE SINVSN.AI , B . 0< 


C 

(-,$***#SURRnUTINE TO GENERATE THE INVERSE OE THE MATRIX A] 


C. 

c 

r, 

c 

c 

r 

r. 

c 

c. 

c. 

c 

c 

C * # # 

c 

c 

c 


the matrices ai and b are stored in general form 
input matrix is ai 

OUTPUT INVERSE MATRIX IS fi 
N IS THE ORDER OF A I 
n IS THE DETERMINANT OF AI 

L - WORK VECTOR OF LENGTH N 
M - WORK VECTOR OE LENGTH N 

THE STANDARD GAUSS-. IORDAN METHOD IS USED. THE DETERMINANT 
IS ALSO CALCULATED. A DETERMINANT OF ZERO INDICATES THAT 
THE MATRIX IS SINGULAR. 

**WITH MODIFICATIONS TO INPUT MATRIX IN VECTOR FORMAT 


DIMENSION A I %400< , R2400< , A2400< , L t 2 DC , M2 2 D< 


C 

C 

c 

c IF A DOUBLE PRECISION VERSION OF THIS ROUTINE IS DESIRED, THE 

C C IN COLUMN 1 SHOULD BE REMOVED FROM THE DOUBLE PRECISION 

C STATEMENT WHICH FOLLOWS. 

C 

C • DOUBLE PRECISION A , 0 , B I GA , HOLD 
C 

C THE C MUST ALSO RE REMOVED FROM DOUBLE PRECISION STATEMENTS 

C APPEARING IN OTHER ROUTINES USED IN CONJUNCTION WITH THIS 


C ROUTINE. 

C 

c THE DOUBLE PRECISION VERSION OF THIS SUBROUTINE MUST ALSO 

C CONTAIN DOUBLE PRECISION FORTRAN FUNCTIONS. ABS IN STATEMENT 

C ID MUST RE CHANGED TO DABS. 

C STORAGE OF AI ELEMENT IN A 

C 


KK#N*N 
DO 5 J # 1 , K K 
5 A*J<#AI£J< 


C 

C 

c 

c SEARCH FOR LARGEST ELEMENT 

C 

D#1 .D 
NK#-N 

DO 80 K# 1 , N 

NK0NKEN 

LXK<#K 

M3tK<#K 

K K # N K E K 

BIGA#A3?KK< 

DO 20 J#K,N 
IZ#N**J-1< 

DO 20 I #K » N 


MOP 1? 630 
MOP], 2 640 
MOP 12 650 
MOP 1 2 660 
MOP] 2670 
MOP 1 2 680 
MOP 12690 
MOP] 2700 
MOP 12710 
MOP 1 2 720 
M 0 P 1 2 7 3 0 
MOP 12740 
MOP12750 
MOP 1 27 60 
MOP12770 
MOP 12780 
MOP 12 7 90 
MOP 1 2 BOO 
MOP 128 10 
MOP 12 820 
MOP 12 8 30 
, MOP 1 2 840 
MOP 12 8 50 
MOP 12860 
MOP 128 70 
MOP 12 880 
MOP 12 890 
MOP 1 2 900 
MOP 129 10 
MOP12920 
MOP 12930 
MOP 12940 
MOP 12950 
MOP 1 2 9 60 
MOP 12970 
MOP 1 2 980 
MOP 12990 
MOP 1 3 000 
MOP 13010 
MOP 1 3020 
MOP 1 3030 
MOP 1 3040 
. MOP 13050 
MOP 130 60 
MOP 13070 
MOP 1 3080 
MOP 13090 
MOP 1 3100 
MOP 13110 
M0P13120 
MOP 13130 
MOP 1 3 140 
MOP 1 3150 
MOP 13160 
MOP 13170 
MOP 13180 
MOP 13190 
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uaizfii 

M0P13200 

10 

IFSARS*BIGA<-ABS*ASI J«< 15,20,20 

MOP 13210 

1 s 

R I GA#A$ I J< 

M0P13220 


L*K<#I 

MOP 1 32 30 


M?K<# J 

MOP 1 3240 

20 

CONTINUE 

MOP 1 32 50 

c 


MOP 1 32 60 

c 

INTERCHANGE ROWS ' 

M0P13270 

c 


MOP132RO 


J#LSK< 

MOP 1 3290 


IFS!J-K<35,35, 25 

MOP 1 3300 

25 

K I #K — N 

MOP 1 33 1 0 


no 30 I # 1 , N 

MOP 13320 


K I I f.N 

MOP 1 3330 


HOLD#-AXKI< 

M0P13340 


J I # K I — K fi J 

MOP 1 33 50 


AXK I <#A3!J I < 

M0P13360 

30 

A*JK#HOLD 

M0P13370 

c 


M0P133R0 

C 

INTERCHANGE COLUMNS 

MOP 1 3390 

C 

' 

MOP 1 3400 

35 

I#MXK< 

MOP 1 341 0 


IF3!I-K<45,45* 3ft 

MOP 1 3420 

3fl 

JP#N**I-1< 

MOP 1 3430" 


00 40 J#1,N 

MOP 1 3440 


JK#NKEJ 

MOP 13450 


JI#JPf.J 

MOP 13460 


HOLD#-ASJK< 

MOP 1 3470 


A*JK<#A$J I< 

MOP 13480 

40 

AXJI<#HOLD 

MOP 1 3490 

C 


MOP135O0 

C 

DIVIDE COLUMN RV MINUS PIVOT (VALUE OF PIVOT ELEMENT IS 

MOP 1 3510 

C 

CONTAINED IN RIGA) 

MOP 13520 

C 


MOP 135 30 

45 

IF*BIGA< 4ft » 46 , 4ft 

M0P13540 

46 

n#o.o 

MO PI 35 50 


GO TO 150 

MOP 1 35 60 

4ft 

DO 55 I#1»N 

MOP 1 3570 


IF*I-K<50,55., 50 

MOP 13580 

50 

I K # N K f> I 

MOP 1 3590 


AXIK<#A«IK</^-RIGA< 

MOP 1 3 600 

55 

CONTINUE . 

MOP 13610 

r. 


M0P13620 

c 

REDUCE MATRIX 

MOP 1 3 630 

c 


MOP 1 3 640 


DO 65 I # 1 , N 

M0P13650 


IK#NK£I 

MOP 1 3 660 


HOLD#A*IK< 

M0P13670 


I J # I — N 

M0P13680 


DO 65 J#1,N 

M0P13690 


I J# I J EN 

MOP 1 3700 


IF35I-K <60,65, 60 

MOP 13710 

60 

IFSJ-K<62,65,62 

MOP 1 3720 

62 

KJ#I J-IEK 

M0P13730 


A*I J<#HOLD*AXKJ<CA*I J< 

MOP 1 3740 

65 

CONTINUE 

MOP 13750 

C 


M0P13760 
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r. divide row by pivot 

r. 

K J#K-N 

nn 75 j # 1 ,n 

KJ#KJfN 

TFXJ-K<70,75,70 
70 A2K J<* AXKJ</RIGA 
75 CONTINUE 
C 

f. PRODUCT OF PIVOTS 

C 

n#n*Rir,A 

c 

c. RFPI.ACE PIVOT by RECIPROCAL 

C 

A2KK<#1 .O/RIGA 
RO CONTINUE 

f, 

C FINAL ROW AMD COLUMN INTERCHANGE 

C 

K#N 

100 K#*K-1< 

I F2!K< 150,150,1 05 
105 I#L*K< 

IFXI-K <120,1 20, 10R 
108 J0#N*SK-1< 

JP#M*£I-1< 

no no j#i,n 
JK# J0 6J 
HOL0#A* JK< 

J I # >) R f. >1 
A2JK<#-AS!JI< , 

110 A*JI<#HOL0 
120 J«MS!K< 

IF3!J-K< 100,100,125 
125 K I #K— N 

on 130 I#1,N 
K I UK I r,N 
HOLD#AXK I< 

.1 1 #K I -K £J 
A*KI<#-A*J1< 

130 A 20) I <#HOL 0 
CO TO 100 

150 L L # 0 
C . . 

KK#N#N 

DO 151 J # 1 , K K 

151 RSJ<#AS!J< 

RETURN 

C 

END 


MOP 13770 
MOP 1. 3780 
MOP 137 90 
MOP 13 800 
MOP 138 10 
MOP 13820 
MOP13830 
MOP 1 3840 
M0P1 3850 
MOP 138 50 
MOP 138 70 
MOP 13 880 
MOP 1 3890 
MOP 1 3900 
MOP 1 39 10 
M0P1 3920 
MOP 1 39 30 
MOP 1 3940 
M0P13950 
MOP 1.3 9 50 
MOP13970 
MOP13980 
MOP 13990 
MOP 14000 
MOP 1401.0 
MOP 1 402 0 
MOP 140 30 
MOP], 4040 
MOP 14050 
MOP] 4050 
MOP 14070 
MOP 14080 
MOP 1 4090 
MOP 1.4 100 
MOP 1.41 ] 0 
MOP 1.43.20 
MOP 1.41 30 
MQP14140 
MOP 14150 
MOP 141. 60 
MOP 1 41 70 
MOP 14180 
MOP 1 4190 
MOP 14200 
MOP 142 10 
MOP 1.4220 
MOP 142 30 
MOP 1 4240 
MOP 142 50 
MOP 14? 60 
MOP 142 70 
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peegtomg page blank not filmed 


APPENDIX D 

INPUT OUTPUT OPERATIONS PACKAGE 

D. 0 Introduction 

This appendix contains listings for the following programs. 


Page 

SUBROUTINE IMXPNP(M«NA.NB»MENTRY) 102 
SUBROUTINE IRAN0R<ND.IA«IB.IC.I0.IE»IF,IG.NENTRY> 103 
SUBROUTINE MXRNP (VA*NA»NB*NENTRY) 104 
SUBROUTINE NAMRNP <M,NA,NB*NENTRY> 105 
SUBROUTINE RANDP (NENTRY) 106 
SUBROUTINE RANOPO (NO «DA «DB *DC *DD * DE *DF »DG» NENTRY) 106 
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SUBROUTINE IMXRNP ! M,NA,NB,NENTRY ) 

SUBROUTINE READS, PRINTS AND STORES INTEGER NA*NB MATRIX 
MATRIX IS STORED IN GENERAL FORM COLUMN BY COLUMN 

DIMENSION M ( 1 ) 

1000 FORMAT! 7110) 

1002 FORMAT! 7115) 

1003 FORMAT! 7F10.4) 

READ IN NA BY NB MATRIX ROW-WISE AND STORE INTO 1 DIMENSION 
VECTOR COLUMN-WISE. 

GO TO! 1,1,2, 4, 2) , NENTR Y 

1 J=NA*NB-NA+1 
DO 15 I - 1 , NA 

READ 1000, ( M< K) ,K=I , J,NA ) 

15 J=J+1 

GO TO! 2, 3, 3, 2), NENTR Y 

PRINT NA BY NB MATRIX ROW-WISE 

2 CONTINUE 
JJ=NA*NB— NA+1 
DO 11 11=1, NA 

IF! NENTRY— 5 ) 12, 10, 12 

10 PUNCH 1003, !M( L > , L= I I , J J ,NA ) 

GO TO 11 

12 PRINT 1002, (M(L),L=II,JJ,NA) 

11 JJ=JJ+1 

C 

3 RETURN 
C 

C READ AND PRINT HEADING CARD BEFORE READING AND PRINTING MATRIX 

4 CALL R ANDP! 4 ) 

GO TO 1 

C 

END 


I0P1Q000 
I OP1O01O 
I 0P10020 
I OP 10030 
I OP 10040 
I OP 10050 
I OP 100 60 
I OP 10070 
I OP 1 0080 
I OP 10090 
I0P10100 
IOPIOHO 
I0P10120 
I0P10130 
I OP 10 140 
I OP 1 01 50 
I0P10160 
I OP 10170 
I OP 10 1 BO 
I 0P10190 
I 0P10200 
1 OP 102 10 
I OP 1 0220 
I OP 102 30 
I 0P10240 
I 0P10250 
IOP10260 
I0P10270 
I0P102B0 
10P10290 
I OP 1 0300 
I0P10310 
I0P10320 
I 0P10330 
I OP 10340 
I OP 10350 
I OP 10360 
I0P10370 
I0P10380 
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SUBROUTINE IRANDPI NO , I A , I B, I C , I 0, I E , I F , I G , NENTR Y ) 

I OP 10390 



I 0P10400 


SUBROUTINE TO READ AND PRINT INTEGER DATA 

I OP 104 1 0 



I OP10420 


DIMENSION I V( 7 ) 

I0P10430 



I OP 1 0440 

1000 

FORMAT ( 7110) 

I OP 1 0450 

1010 

FORMAT ( 7115) 

I OP 104 60 



I0P10470 


GO T0{ 1,2, 1,4) ,NENTRY 

I OP 1 0480 

1 

READ 1000, I A , I B , IC, ID, IE, IF, IG 

I OP 10490 


GO T0( 2,2,3,21 ,NENTRY 

I OP 10500 

2 

I V( 1 ) = I A 

I OP 105 1 0 


I V( 2 ) = I B 

I0P10520 


IV( 3 ) = I C 

I OP 105 30 


I V( 4 ) = I D 

I OP 1 0540 


IV< 5) = IE 

10P10550 


I V( 6 ) = I F 

I 0P10560 


IV( 7)=IG 

I0P10570 


PRINT 1010, ( I V( I ) , 1 = 1, ND) 

I OP 1 0580 

3 

RETURN 

I OP 1 0590 



I0P10600 

4 

CALL RANDPI4) 

I OP 1 06 1 0 


GO TO 1 

I 0P10620 



I OP 10630 


END 

I OP 1 0 640 
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SUBROUTINE MXRNPI VA ,NA ,NB ,NENTRY ) 


I OP 10650 


SUBROUTINE READS AND/OR PRINTS 
GENERAL FORM COLUMN BY COLUMN 


THE NA*NB MATRIX VA WHICH IS 


I 0P10660 
ST0REDI0P10670 
I OP 1 0 680 


DIMENSION VA( 1 ) 

1000 FORMAT ( 7E10.0) 

1002 FORMAT ( 7F15.6) 

1003 FORMAT! 7110) 

GO T0( 1,1,2, 4,2), NENTRY 

READ IN NA BY NB MATRIX ROW-WISE AND STORE INTO 1 DIMENSION 
VECTOR COLUMN-WISE. 

1 J=NA*NB-NA+1 
DO 15 1=1, NA 

READ 1000, (VA(K),K=I,J,NA) 

15 J=J+1 

GO T0( 2,3,3,2),NENTRY 

C 

2 CONTINUE 

C PRINT NA BY NB MATRIX ROW-WISE 

JJ=NA*NB-NA+1 
DO 11 11=1, NA 
I F ( NENTRY-5) 12, 10, 12 

10 PUNCH 1003, (VAIL), L=II,JJ,NA) 

GO TO 11 

12 PRINT 1002, (VA( L)»L=II,JJ,NA) 

11 JJ=JJ+1 
RETURN 

C 

3 RETURN 
C 

C READ AND PRINT HEADING CARD BEFORE READING AND PRINTING MATRIX 

4 CALL RANDPI4) 

GO TO 1 

C 

END 


I OP 1 0690 
I OP 10700 
I OP 107 10 
I OP 10720 
I OP 1 07 30 
I OP 1 0740 
I0P10750 
I OP 1 07 60 
I0P10770 
I0P10780 
I0P10790 
I OP 1 0800 
I0P10810 
I0P10820 
I OP 10830 
I OP 10840 
I0P10850 
I OP 1 08 60 
I0P10870 
I 0P10880 
I0P10890 
I 0P10900 
I OP 1 09 1 0 
I 0P10920 
I0P10930 
I OP 10940 
I OP 10950 
I OP 1 0960 
I OP 1 0970 
I OP 10980 
I OP 1 0990 
I0P11000 
I0P11010 
I0P11020 
I OP 1 1030 
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SUBROUTINE NAMFUVP! M,NA,NB,NENTRY) 

SUBROUTINE READS, PRINTS AND STORES INTEGER NA*NB MATRIX 
OF FOUR CHARACTER NAMES 

MATRIX IS STORED IN GENERAL FORM COLUMN BY COLUMN 
DIMENSION M( 1 ) 

1000 FORMAT ( IX , A4, IX , A4 , 1X,A4,1X»A4, 1 X » A4, IX * A4, IX , A4, IX * A4* IX, AA, 
1 1X,A4,1X, A4, IX , A4, IX , A4, IX , A4) 

1002 FORMAT! 1 1 X , A4 , 1 1 X , A4 , 1 IX , A4, 1 IX , A4, 1 1 X , A4 , 1 IX , A4, 1 1 X , A4 ) 

GO TO! 1, 1 ,2,4) ,NENTRY 

READ IN NA BY NB MATRIX ROW-WISE AND STORE INTO 1 DIMENSION 
VECTOR COLUMN-WISE. 

1 J=NA*NB-NA+1 
DO 15 1=1, NA 

READ 1000, (M(K),K=I,J,NA) 

15 J=J+1 

GO TO! 2, 3, 3, 2 ) ,NENTRY 

PRINT NA BY NB MATRIX ROW-WISE 

2 CONTINUE 
JJ=NA*NB-NA+1 
DO 11 11=1, NA 

PRINT 1002,! M< L ) ,L=I I , JJ,NA) 

11 JJ=JJ+1 

C 

3 RETURN 
C 

C READ IN HEADING CARD BEFORE READING AND PRINTIND M 

4 CALL R ANDP! 4 ) 

GO TO 1 

C 

END 


I OP 1 1 040 
I OP 1 1050 
I OP 1 1060 
I OP 1 1070 
I OP 1 1 080 
I OP 1 1090 
I0P11100 
IOPllllO 
I0P11 120 
I OP 1 1 1 30 
I OP 1 1 140 
I OP 1 1 1 50 
I OP 1 1 1 60 
I0P11170 
I0P11180 
I OP 1 1 1 90 
I OP 1 1 200 
I0P11210 
I OP 1 1 220 
I OP 1 12 30 
I OP 1 1240 
IOP11250 
I OP 1 1 2 60 
I0P11270 
IOP11280 
I OP 1 1 290 
I OP 1 1 300 
I OP 1 1 3 10 
I OP 1 1 320 
I OP 1 1 330 
I DP 1 1 340 
I OP 1 1350 
I OP 1 1 360 
I0P11370 
I OP 1 1 380 
I OP 1 1 390 
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SUBROUTINE RAN[)P( NENTRY ) 

I OP 1 1 400 

c 


I OP 1 1 A 1 0 

c 

SUBROUTINE TO READ AND PRINT HEADING CARDS 

IOP11420 

c 


I OP 1 1430 


DIMENSION FNAMEI fl) 

I OP 1 1 440 


DOUBLE PRECISION FNAME 

I OP 1 1450 

c 


I OP 11460 

1000 

FORMAT ( BA8) 

I0P11470 

1001 

FORMAT ( 1 H 1 ) 

I OP 1 1480 

1010 

FORMAT ( 2X,A8,2X,A8,2X,A8,2X,A8,2X, A8,2X,A8,2X,A8) 

I0P1 1490 

1020 

FORMAT ( 7X,A8,7X,A8,7X,A8,7X,A8,7X,A8,7X,A8,7X,A8) 

I OP 1 1500 

c 


I0P11510 


GO T0( 2000,2010,2020,2030) , NENTRY 

I OP 1 1520 

c 


I OP 1 15 30 

2000 

PRINT 1001 

I OP 1 1 540 

2010 

READ 1000, (FNAMEI I), 1 = 1, 8) 

I OP 1 1550 


PRINT 1000, (FNAME( I), 1 = 1, 8) 

I OP 1 1 560 


RETURN 

I0P11570 

c 


I OP 1 1 580 

2020 

PRINT 1001 

I OP 1 1590 

2030 

READ 1010, ( FNAMEI I ), 1=1,7) 

I OP 1 1 600 


PRINT 1020, ( FNAMEI I ), 1 = 1,7 ) 

I0P1161G 


RETURN 

I OP 1 1 620 

C 


I OP 1 1630 


END 

I OP 11 640 



SUBROUTINE RANDPDI ND , D A , DB , DC , DD , DE , OF , DG , NENTR Y ) 

I OP 1 1 650 

C 


I OP 1 1 660 

C 

SUBROUTINE TO READ AND PRINT FLOATING POINT DATA 

I OP 1 1 670 

C 


I OP 1 1 680 


DIMENSION DVI 7 ) 

I OP 1 1 690 

C 


I OP 1 1 700 

1000 

FORMAT! 7E10.0) 

I0P1 1710 

1010 

FORMAT! 7F15.6) 

I OP 1 1720 

C 


I OP 1 1730 


GO TO! 1,2, 1,4), NENTRY 

I OP 1 1740 

C 


I0P11750 

1 

READ 1Q00»DA,DB,DC»DD,DE,DF,DG 

I OP 1 1 7 60 


GO TO! 2, 2,3,2) , NENTRY 

I OP 1 1770 

2 

DVI 1 )=DA 

I OP 1 1 780 


DVI 2 ) =DB 

I OP 1 1790 


DVI 3 ) =DC 

I OP 1 1 800 


DVI 4)=DD 

I OP 1 18 10 


DVI 5 ) =DE 

I OP 1 1 820 


DVI 6 ) =DF 

I0P11830 


DVI 7 ) =DG 

I OP 1 1 840 


PRINT 1010,1 DVI I ) , I = 1 ,ND ) 

I OP 1 1850 

3 

RETURN 

I OP 1 1 860 

C 


I0P11870 

4 

CALL RANDPI4) 

I OP 1 1880 


GO TO 1 

I0P11890 

C 


I OP 1 1 900 


END 
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